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Summary 


Some essential parts of the spray drying process are treated, viz. 


1. The evaporation of liquid drops, both during the first period after formation when they are decelerated 
by air resistance, and during the second period when they fall at a constant settling rate. 
2. The humidity and temperature of the drying air during the process with reference tu the humidity 


chart. 


The transfer units-concept is introduced and a rational evaluation of the spray drying process is suggested. 


Zusammenfassung —Es werden einige essentielle Teile der Zerstdubungstrocknung besprochen und zwar: 
1. Die Verdampfung von Tropfen, sowohl wahrend der ersten Periode nach ihrer Formierung, wenn sie 


durch Luftwiderstand gehemmt werden, als in der zweiten Periode wenn sie mit konstanter Endgeschwindig- 


keit fallen. 


Feuchtigkeitsdiagramm. 
Der Begriff der “transfer units” 
Zerstaubungstrocknung wird vorgeschlagen. 


1. INTRODUCTION 


The object of spray drying is to recover a solid mat- 
erial, dissolved or suspended in a liquid, from this 
liquid in the form of a dry powder, without damaging 
it thermically during the operation. This entails that 
the operation in general is of short duration (a few 
seconds) and that it must be carried out at moderate 
temperatures. 

The process as such consists in the contacting of 
a dry, hot gas with drops produced by atomization 
from the solution or suspension in a thermally insu- 
lated space, the drying chamber. During the contacting 
the liquid of the drops evaporates, while the gas be- 
comes moister and colder. The mixture leaving the 
drying chamber consists of a powder and a gas, which 
are separated by means of a cyclone and/or filters. 
The contacting process of drops and gas is often carried 
out under concurrent conditions. 

There are many aspects to the spray drying process. 
The following subjects will not be treated in the pre- 
sent article: the air heater; the drying chamber; the 
droplet-formation; the dropsize-distribution; a com- 
parison with other drying processes; economical 
considerations. 
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2. Feuchtigkeitsgrad und Temperatur der Trocknungsluft wahrend des Prozesses in bezug auf das 


wird eingefiihrt und eine rationelle Methode zur Beurteilung der 


We shall, however, contribute some theoretical 
views on some of the essential parts of the spray drying 
process, 

These refer to the evaporation process of the drops 
in moving air. Primarily it is shown that there is a 
marked correlation between the heat and material 
transfer in the gas phase. A calculation can be made 
of the evaporation of the drop, both during the first 
period after formation when the drop is decelerated 
and during the second period when the drop falls at 
a coustant rate. 

Further, the change in condition of the gas during 
the process will be discussed with reference to the hum- 
idity chart. As in practice the liquid to be evaporated 
is very often water and the gas used is air several of the 
formulas given and of the conclusions arrived at apply 
to this combination. 


2. THe Evaporation oF A Drop oF Liguip IN AIR 
The simplest case is that of a spherical drop of pure 
liquid of diameter D and temperature 7’ in stationary 
air. The vapour tension of the liquid is P(7'), the 
partial vapour tension in the air, p, being smaller 
than P(T). 
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As the vapour concentration at the surface of the 
sphere is higher than in the air atmosphere, there is 
molecular diffusion of vapour molecules outwards *. 
Lanemurr [1] calculated the total diffusion current 
for this case, viz. 

N —2aDD Ip kmol 


"RT sec 


(1) 


A more general problem is that of a drop which 
travels through the air at a certain velocity r. Vapour 
removal will then be promoted, because the vapour 


molecules only have to diffuse through a thin film 
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In McApams [3] these data are correlated as the rela- 
tion between the dimensionless numbers of Nusse.t 
and REYNOLDs: 


, AD eDe 
N : 
Nu ky / u ( 3) 
FROSSLINGs evaporation equation 
, Ip 70-5 y70-33) kmol 
2a DD (276 Np, Ng 2 
N aDD, RT l 16 Nae Nse | (2) 


can be replaced by the corresponding equation for 
heat transfer: 


dw J 








: ° ’ a ove NOS yO 33 
of air round the drop. dg AAD, AT | + O26 Se hr? 
Naw - i 
Nussevts number can be derived from 
_ (4). When at the same time the Pr-value 
my §.. for air is substituted the result is: 
« é 7 
Roo?* | , AD ’ an 10'S 2 
° to eet Nyw ky sd O25 Nae. (5) 
10a} Po) ut ; 
* a * eee ee : 
:= Eq. (5) is in good agreement with 
e . 0.5 ’ ; : ’ 
g Nn 29.05 Nig Witiiams’ graph in McApams’ book men- 
M tioned above (Fig. 1). This proves the 
reliability of FrOssLines equation. 
For the case of evaporating drops of 
water in atmospheric air eq. (2) gives on 
3 substitution of the following values: 
5 : , 
| . D,— 30.10 8m*® see, Ny = 0-60, 
10 10? 10? 10° 10° Nags 0G 1-28 kg m®, M IS 
Fig. 1. Exp. data for heat transfer from air to spheres. Source: and 
McApams’ “Heat Transmission.”’ Theoretical line, derived from 
FROSSLINGs data. Ip a in 1H CG 1H 
RT . ‘M I4-1 : 
(6) 
According to FrOssLING [2] the following empirical — vith Ht in ke H,O 
: s F ‘ a kg dry air 
equation is valid for the evaporation of drops in mov- : 
ing air \ Ia BO. 10 4° D 1H 
— me ra) 
, Ip =~ 05 yy0-33) kmol 
» | -276 N, » ne wet ~' kmol 
N= 32DD, pg +088 Me): 1 + 0-276 Na” 0-843 ° 
sec 
Froésstinc used nitrobenzene, aniline, water and 


naphthalene (solid). The diameter of the drops was 
(-2-2 mm, the relative air velocity 0 to 7 m see. 

In eq. (2) Np, represents the influence of the 
velocity on the condition of flow round the sphere, 
hence on the thickness of the film. Scumipts number 
gives the influence of diffusion in the same way as 
does PRANDTLs number in the analogous problem of 
heat transfer. 

Hence an additional justification of eq. (2) for the 
evaporation of a drop moving in air may be found 
by comparison with data on heat transfer to spheres. 


* The case of natural convection was treated by P., Meyer; 
Trans. Inst. Chem. Eng. 1937 15 127. 


By multiplication by the molecular weight of water, 
M Is: 


L =2-4.10°*DAH 1 4+ 0-23 Np,” kg/see. (7) 


Fractional evaporation follows from this by division 


by the weight of a drop, W - D3, 1000. 


d x, ; - 4H oe ys). , 
de 46.10 pe 1 + 0-23 Np,” sec”!. (8) 
Though in a spray drying plant the movement of 
the drops is, admittedly, highly complicated and dif- 
ficult to express in formulas, the braking of the high 


+ H = 0-001 about corresponds to 1-2mm Hg partial 
water vapour pressure. 
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initial velocity and then the free fall in relation to the — the water of its own accord assuming the wet bulb 
moving air are two of its essential elements. Therefore, temperature of 43°C. This is further explained in 


while bearing in mind this schematic picture we will — section 7. 


“Re, vin m/sec FRACTIONAL EVAPORATION 


#4.6.107 Att (i¢0.23 Ne? ) 


aN 
2) te *4 42 an [Sts 
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Fig. 2. Constant settling velocity and Re-Number of free Fig. 3. Theoretical equations for fractional evaporation 
falling water drop in air. Odrop 1000 kg m®, air of water drops in air. 1) Free falling drop (per metre) 
Lkg/m?*, sir 2-1. 1075 kg/mseec, ¢ 83°C. 1H — 0-06; 2) Brake-evaporation (total) 4H 0-06. 


now deal with the evaporation of a 
drop during these two stages. 


Tue Free Fatiine Dror 
or Liguip 
We apply the equation found above 
(8) to a falling drop of water that 
has reached its final velocity v in 
air. This velocity greatly depends 
on its diameter (see Fig. 2). 

The fractional evaporation per 
metre length of fall (relative to the 
air) follows by division by the rela- 
tive velocity = v: 
dx . 


dl - 


Fig. 4. Exper. values of the drag coefficient C, as a function of the REYNoLDs’ 
number Np,, for spheres in a fluid. 


Lo 


1H r0-5 

7 23 2 

l » Dt l 3 Np, 

To apply eq. (9) at a given diameter D one has to Calculation of the constant settling speed of a falling sphere 
find the corresponding values of v and Np, and then (compare Perry [4}) 

; » . i € 5} ’ i 4 , i . . . . . 

substitute them in (9) at a certain value of 4H, which — Ppo resistance to the motion of a sphere in a fluid can be ex- 
follows from the temperature of the drop and the pressed as a dimensionless drag coefficient C,, which is plotted 
humidity of the air. In Fig. 3 the fractional evapora- against the ReyNnotps number (see Fig. 4). The drag coeffi- 
tion for drop diameters between | and 0-1 mm is __ cient C, is defined by 
plotted at 1H = 0-06. This 4H is reached, for in- 


stance, in the case of water drops in dry air of 193° C, (10) 
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For falling spheres also 


F 9. Pins 0). (11) 


For constant settling velocity, (10) is equal to (11), so 


49 D (os 0) 











Chemical 
Engineering Science 


of account, the velocity is braked by the fluid. If the 
final speed is given, calculation of the distance covered 
as well as of the braking time is required. 


In retarded motion the kinetic energy of the sphere is 


partly consumed to overcome the resistance. This is expressed 


30a by the equation: 
or 
. 4q LF o(0.— 0) d l . > dv > r 
Cy Np . SSS, 12 me" Fe or m F. 13 
dN Re 34 (te) dr \2 ' dt (13) 
{8 Now 
! Re 
Ie ee © a 
fk Cd. » 2G" 4 1» and m 6 IP .o,- 
aan Therefore : 
10 Pe 
1 dv l 
: ~ DP 4 -.C,.05.0°. BP. l4 
~ 3 CL dr 4 d Cg ( ) 
: ‘ . 
= Now introduce a dimensionless velocity and time: 
, Div T 
| Np, and (a) - 
0" of 4o, BP 
aNae 68 3 0G -— 
2 
argh frac Then the eq. (14) reduces to: 
d iN VOL, 
10 de a ; vi ' (15) l 
ds Ke 1952 
is The time 11 which elapses when the drop is slowed 
down from the velocity v, to vy, i.e. for a change of Re- 
number from (Np), to (Ny.),, is: 
ee eee eee : (Npe); 
Fig. 5. Moving sphere, braked by a fluid Ir 4o, BP Ja e 40, BP | d Nr (16) 
. ad . . ° _ . ‘2 
Graphical determination of reduced time @ 30% *,. 3% Cy. Ni, 
("pede 
ior , 2? . . . y 7 wT é ’ . - ‘ ’ d N le . 
In Fig. 4, Cy Np, is plotted as a function of Nye. Fora given We call | do | Kt the reduced time @. It is plotted 
sphere, by eq. (12) we can calculate Cy Np,. From Fig. 4 we ¢ J Ca. Ni, 


find NV, — e Deu, so v can be determined. 


4. Tue Brakep Drop or Liguip 
Before the drops in a spray drying plant have reached 
their final velocity (discussed above) they have often 
had a much higher initial velocity. This happens, 
for instance with the much-used rotary dise atomizer, 
where the drops leave the dise at a very high horizontal 
velocity (say 150m sec). The drops are very soon 
slowed down by the air, but this is accompanied by 


considerable evaporation. 


(a) Path of the drop 
A precise method is now given for calculating the brak- 
ing distance and braking time for any drop in any 
medium in any velocity interval. 

A sphere is shot into a gas or liquid at a certain 
initial velocity. If the influence of gravity is left out 


in Fig. 5. 


The distance .1/ covered during this time Ar is given by: 
Al= fv.dr. 

1 
We again introduce Np, and 0. Then: 


2 (Ype), 
4,@¢ fy 4 CL “  ANey a 
Ti ) Nr, .dO » » ‘ —— o 
3 d 0G | hed 3 i % | Ca. Nay (l ‘) 
l (Np, ls 


, P as * aANn , 
We call | Np, dO | 7 Ne the reduced distance o. It 


e d+ ¥Re 
is plotted in Fig. 6. 
Special cases: 


(a) Case of Newtons law, C, > 0-44, obtaining for 700 
Np, — 200000, 


It is easy to calculate from eq. (16) and (17) that for this case 


Ar = 302 “4 p : 18 
: (15) 
CG 3 %, | 
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Al = 3-02 en Din { 2 |. 
V2 


Cg 


(b) Case of Stokes’ law, Cy = 24/Np,, for Np, < 14. 


For this case we find in a similar way from eq. (16) and (17): 
Ar = 0-0555 


cL [In vr, —In vy) (20) 
0g 


Al = 0-0555 (21) 
Example of the use of graphs 5 and 6 to determine the braking 


time and braking distance 


Let us suppore that a drop of water with an initial velocity 
of 67 m/sec is braked by the air to a final velocity of 5 m see. 
The diameter of the drop is 100 u. The densities of water and 
air are related as 0, 0, — 800. The kinematic viscosity of 
air is y = 15.1075 m? see. 

We first calculate the initial and final values of 
446 


and find Np, and N 33. 


In graph Fig.5 we find 10 — 0-031 — 0-005 — 0-026, 
Therefore 
4 PP 


‘ O-O1S85 sec . 
3 


Ir 


In graph Fig. 6 we find Jo 30-3 = 3-4. 
Therefore 
: °L lo = 0-36 m. 
CG 

LarreLe and SuHerverp [5] give extensive calculations for 
four cases: 

1. Two-dimensional motion in a gravitational field. 

2. One-dimensional motion in absence of a gravitational 
field. 

3. One-dimensional motion in the presence of a gravi- 
tational field. 

4. Two-dimensional motion in a centrifugal field. 

Of these, No. 2 partially covers the above-described case : 
the braking time integral (16) Fig. 5 is found, the braking di- 
stance integral (17) Fig. 6, is not, however. 


(b) The evaporation of the drop 


When a drop is shot into a gas at a high initial speed 
considerable evaporation takes place before this 
movement is entirely braked. It was discussed above 
how this braked movement exactly proceeds. With 
the aid of the equations developed and with Frdéss- 
LINGs evaporation equation it is possible to calculate 
this evaporation during braking. FRrOssLINGs equation 
reads: 


N= 22 DD, 4? (1 + 0-276 Ng No) Kol 


CRT S sec 


Spray drying 


In this equation the second term represents the extra 
evaporation because of the movement of the drop. 
Take D to be constant. 

Then the extra evaporation during the braking 
period is given by: 


, 1 _ 70 33 . 705 
E,=22DD, 7", .0-276 Ngo { Np dt. (22) 
1 


: "ee 








. - | | 

+ ; 

J 

oe w T* fries ® 
Moving sphere, braked by a fluid. 

determination of reduced distance oa. 





Fig. 6. Graphical 


We again introduce the reduced time 9. Then, using 
15): 
(Npe), 
| * _d Nn 
’ ri5 * 
. ! d* Nie 
(Nye le 


Dd, Ip 
“RT 


. a a ()- 3! 
E. = 0-276. = De nos 23) 


The fractional evaporation is found by dividing 

iP ~AH.%, 80 
RT = 
(Nre), 


4-42 AH Nz? [ 


77 o 
pel 
6 ! M’ 


and as 
- _ ? 24) 
ned a+ Re 

It may be noted that the fractional evaporation is 
independent of the drop diameter and that, except for 
some material constants, it only depends on the initial 
and the final Re-number of the drop. For high initial 
Re-numbers the braking evaporation only depends on 
the final Re-numbers, as shown in Fig. 7. If the con- 
stant settling speed of the liquid sphere is taken as 
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an obvious final velocity of the braking by air, the 
final Re-number depends on the diameter of the sphere, 
as shown in Fig. 2. In Fig. 3 the total evaporation for 
this case is shown, calculated with eq. (24) and Fig. 7. 
It was assumed that the drop diameter is constant —so 
this applies when the fractional evaporation is small, 
177 is small. 


If, however, the diameter may not be assumed to 


eq. because 


be constant. the braked movement is different and the 
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Fig. 7. Moving sphere, braked by a gas. Graphical 


determination of the extra evaporation 


extra evaporation must be ascertained step by step; 
with decreasing diameter it will become smaller than 
indicated by eq. (24). 

This is evident, as eq. (24) indicates the fractional! 
evaporation at every moment as referred to the drop 
mass if the diameter should not decrease after that 
moment. In reality the diameter does decrease, so the 
velocity decreases more rapidly than is indicated in 
eq.(14). So the time of evaporation becomes shorter 
than when the diameter is constant. Moreover, a 
smaller area is then available for evaporation than 
when the diameter is constant. Since the two factors, 
time and area, are smaller, the total evaporation as 
referred to the original weight of the drop, will be 


smaller than when the diameter is constant. 
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5. Some CONCLUSIONS FROM THE EQUATIONS FOUND 
FOR THE FREE FaLuLInG AND Brakep Drop or Liguip 


We shall draw some conclusions from the equations, 
as these apply in the first place to the evaporation of 
pure water, 

Fig. 3 shows the fractional evaporation as a fune- 
tion of the drop diameter. The fractional evaporation 
per metre of height of fall greatly decreases with in- 
creasing diameter. 

1. This means that even a comparatively slight 
inhomogeneity of drop sizes immediately results in 
the smaller drops drying too rapidly (with the chance 
of scorching) or the larger drops too slowly (with the 
chance of too great water content) 

2. Another conclusion drawn from this curve is 
that moderate reduction of the drop size causes a great 
acceleration of the fractional evaporation, for which 
a much smaller height of fall is required. When, for 
instance, the required drop diameter is reduced from 
0-9 to 4mm the fractional evaporation rate is in- 
creased from 3-6 to 24%, so that the height of fall can 
also be decreased far below its original value. 

3. It also follows from eq. (9), which describes the 
curve concerned, that for constant diameter drops the 
fractional evaporation is approximately in inverse 
ratio to the constant settling velocity when Np, is 
small. Here the special advantage of hollow drops in 
the process becomes evident (moreover, a product 
consisting of hollow drops often has great advantages 
as regards easy solubility and slight bulk density). 
For instance, suppose that only a third of the volume 
of a hollow drop consists of a fluid, two thirds being 
vas. The rate of fall will then be about three times as 
low as in the case of a solid sphere of the same diameter. 
Moreover, evaporation on a weight basis will be equal, 
therefore on a fractional basis three times more rapid 
than in the case of a solid sphere of liquid at the same 
rate of fall. Therefore the “fractional” evaporation 
per metre of height of fall is about 3.3 < 9 times as 
great. 

Fig. 7 shows, that the fractional evaporation dur- 
ing the total braking period is independent of the 
initial velocity (provided the latter is high) but is a 
function of the final ReYNoLps-number only. 

From Fig. 3 the following conclusions can be drawn: 


1. With small drops the “braking evaporation” 
is insignificant in comparison with free falling evap- 
oration, but with larger drops the braking evaporation 
must certainly be taken into account. For instance, in 
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the case of a drop of 0-8 mm diameter the fractional 
evaporation per metre is 4-6% , the total braking evap- 
oration being 13-6%, so the braking evaporation is 
here equivalent to 3m free fall (of course, a certain 
distance is covered in braking evaporation—but this 
may be mainly horizontal, as with a rotary disc 
atomizer). 

2. The fact that the braking evaporation is inde- 
pendent of the initial velocity (provided the latter is 
high as referred to the final velocity) has as its conse- 
quence that there is little point in choosing too rapid 
an initial velocity as, though this considerably in- 
creases the braking distance, it has little or no effect 
on the braking evaporation. Of course, it remains 
possible that the initial velocity should be very high 
for some other reason, for instance to increase the 
capacity of the atomizer or to obtain a suitable drop 


diameter. 


6. Tue Source or THE Evaporation Heat 
oF THE Drops 

In theory, evaporation heat can be derived from three 
sources : 

1. heat content of the drop, 

2. radiation heat from the wall of the drying 
chamber, 

3. sensible heat of the air. 
Rel 
a minor part as its evaporation heat is very high 
(about 2250kJ/kg) as compared with the heat content 
(specific heat is about 4-2 kJ/kg “C). 


In the case of water the first source plays 


Re 2—In a technical plant radiation scarcely makes 
itself felt, as every drop is surrounded by others hav- 
ing about the same surface temperature. In a narrow 
tube this does not apply, however, and allowance can 
be made for this influence. It is not very great; for 
instance, for a 100% radiating wall of 227° C we find 
for water drops of 0-6 mm in air, also of 227° C, that 
the heat transfer is increased by about 4% by radiation 
effects. As the influence of radiation is so small, we 
shall now neglect this effect in our calculations. 


Re 3—In practice, therefore, evaporation heat is 
derived from the sensible heat of the air, which is only 
transferred to the drop when this has obtained a 
sufficiently high difference of temperature from that 
of the surrounding air, by self-cooling. This self- 
cooling automatically adjusts itself in such a way that 
the difference in temperature becomes just sufficiently 
large to provide the necessary evaporation heat by 
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heat flow to the drop, which is the result of molecular 
heat conduction through the air film round the drop. 
This means, therefore, that when air of a certain 
temperature and humidity is placed in contact with 
a drop of liquid, the latter will very rapidly assume a 
certain temperature (wet bulb temperature), which 
generally depends on the liquid chosen, the humidity 
and temperature of the air. 

The air at the surface of the liquid is saturated, so 
here the air point is on the saturation curve. The 
accompanying temperature can be determined as 
follows, see Fig.8. The heat flowing towards the sur- 
face of the drop must be = (liquid molecules diffusing 
(evaporation heat): 


Ak, AHA or 


away) 
1H h 


- 
n“"La ™ 


AhAt 

On the other hand one can write an eq. (26) for the 

adiabatic cooling line as the condition that loss of 

sensible heat of air = evaporation heat of absorbed 
liquid vapour (both as referred to 1 kg dry air): 


C,dt=2.dH or OT =e. (28) 
dt / 
1H 


d ; 99 
and are equal when the “‘psychrometric ratio 


H 
It dt 


h 
Cy. kun 

and, according to Lewis [6], this happens to be the 
case almost exactly with air and water. This means 
that for this combination the wet bulb point is found 
from a certain air point L by passing along the adia- 
batic line until the saturation line is intersected. During 
the whole period of evaporation the wet bulb tem- 
perature remains practically constant, and equal to the 
adiabatic saturation temperature. For other liquid /gas 
combinations this simple method cannot be applied. 
For these it is necessary to determine separately the 
wet bulb point and the adiabatic saturation tem- 
perature, with the aid of eq. (25) and (26) respectively. 
Moreover, upon further wetting of the gas, the wet 
bulb point changes in the direction of the adiabatic 
saturation temperature, with which it eventually 
coincides. 

Reverting to the water/air combination, it may be 
observed that the identity of (25) and (26) applies, 
independent of the value of the evaporation heat /. 
This means that when the latter changes (for instance, 
in the case of a salt recrystallizing, 2 must be reduced 
by the heat of crystallization released, see WILLIAMS 
[7], or, when almost dry material is used 2 must be 
increased by the adsorption heat of the water), the 
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method for finding the wet bulb point remains valid 
all the same: draw a tangent in the point of the 
adiabatic line where the wet bulb temperature is 
required, and intersect this tangent by the vapour 
pressure line. It must be borne in mind, however, 
that an increase of 2 entails a decrease of dH/dt, so 
of the slope of the adiabatic line. This is therefore 


a 
kg aiR 





‘oe 


WBT AST 





oBT 


90% WATER CONTENT 
30 








Fig. 8. A) Scheme of humidity, chart.; B) Scheme of gra 
phical determination of the driving force 4H at 90 and 
30% water content. 


accompanied by a simultaneous shifting of the wet 
bulb point along the vapour tension line, resulting in 
a lower wet bulb temperature. The reverse obtains 
when A becomes smaller. 

Another change generally occurs when solutions 
the 
decreased upon progressive reduction of the humidity. 


and suspensions are dried: vapour tension is 
This (lower) vapour tension line takes over the func- 
tion of the vapour tension line of pure liquid in deter- 
mining the wet bulb point. This of course affects the 
driving force (and hence the rate of evaporation). 
Generally this becomes smaller owing to the reduced 
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vapour tension curve and with increasing evaporation 
heat. Fig. 8 illustrates this for a case in which the 
vapour tension curve is lower and the evaporation 
heat is higher for 30°% water content than for 90% 
water content, 

We emphasize the fact that this determination of 
the wet bulb point only applies to water and air, but 
that generally for other liquid/gas combinations the 
wet bulb point cannot be found with the aid of the 
adiabatic line, nor is it constant during the entire 
period of evaporation. 

As we have already said, the wet bulb temperature 
is determined by molecular heat conduction on the 
one hand and molecular material diffusion (evapora- 
tion) on the other, through the same air film round the 
drop. The two processes are therefore fundamentally 
linked together, as one cannot occur without the other. 
The consequence of this is that we can test Fross- 
LINGs eq. (2) for the evaporation of moving drops of 
liquid by already familiar measurements of heat 

transfer of moving spheres in air. As was further 
explained in section 2, FrOssLines eq. (2) can be 
transformed to 
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Ny u 2\1 25 Nre } . (5) 


Another consequence of linking heat and material 
transfer is that the calculation of the spray cooling 
process can be carried out in a way almost identical to 
that for a spray drying process. Naturally, in spray 
cooling cold air is used, which withdraws heat from 
the drops (sensible heat, crystallization heat, solid- 
ification heat). 

7. Tue Hoemiprry Cuarr 

The evaporation of drops having been discussed in the 
previous sections, we shall now ascertain what changes 
the drying air undergoes during the process. As the 
air changes both as regards temperature and humid. 
ity, the humidity chart (modified MoLurer diagram) 
can obviously be used for the purpose; this gives us 
an air point by plotting against each other the ab- 
solute humidity H (lb H,O per lb dry air) and the 
normal (dry bulb) temperature. See Fig. 9. 

Various quantities can be found in this humidity 
chart, such as dew point temperature, the saturation 
line, the percentage absolute humidity, the percentage 
relative humidity, the adiabatic cooling line and the 
adiabatic saturation temperature. The reader is referred 
to Perry [4]. 

Fig. 10 gives a schematic diagram of the course of 
the evaporation process. The saturation pressure of 
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water is shown by the line 
Sat. If the condition of the 
inflowing air may be re- 
presented by point L, then, 


umd HEAT 


during the evaporation of 
the water, the condition of 
the air changes along a 
known line (the adiabatic 
cooling line) which is given 
by the condition that the 
evaporation heat is supplied 
entirely by the sensible 
heat of the air. This line 
is almost straight and has 
a slope of dH/dt = about 
1/2300° C. This is valid 
for all the adiabatic lines 
of this chart. 

As already shown (see 
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Fig. 9. Humidity chart for low temperature range. The adiabatic cooling lines on this 


section 6) the adiabatic line chart are straight and parallel. From ‘Elements of chemical engineering,’ 2nd ed., 
can also be used (at least by Bapeer and McCane. Copyright, 1936. Courtesy of McGraw-Hill Book Co., New York. 


in the case of the air/water 


system) to find the wet bulb temperature; for this 4 


is given by the point of intersection W of the adiabatic 
line and the line of saturated vapour tension. We see, 
therefore, that in this system the wet bulb temperature 
of pure water does not change in accordance with the 
condition of the air. The driving force, expressed 
either as 41H or as At, decreases, however, as the air 
point shifts along the adiabatic cooling line towards 
the wet bulb point W. The rate at which the air point 
shifts is, at a given weight of evaporated water, in 


inverse ratio to the quantity of drying air. 
Decreasing vapour tensions occur when solu- 
tions or suspensions are dried ; not with pure wa- 
ter. According as the water content of the drop 
decreases, the vapour tension (at the same tem- 
perature) drops. When this dependence is repre- 
sented schematically, as in Fig. 11, we see that 
as the air point rises along the adiabatic line. 
the wet bulb point descends along the adiaba- 
tic line (contrary to what happens with pure 
water, see Fig. 10). Consequently, the wet bulb 
temperature is no longer constant, but rises. 
The relation of the shifts of the air and wet 
bulb points is in direct ratio to the weight 
relation of the currents of liquid and air. 











Fig. 10. Evaporation of water (wet bulb temperature 
constant). 


100% 70% 40% 20% 10% WATER 


CONTENT 





= = 8 oe 4321 
WET BULB TEMP = ORY _ BULB TEMP 
Fig. 11. 
Drying of a wet material (increasing wet bulb temperature). 


The second consequence of the shifting of the wet evaporated. When the point of equilibrium has been 
bulb point is that the driving force (4H or At) de- reached (in an adiabatic process of infinite length), 
creases much more rapidly than when pure water is the moisture content of the drop becomes so small 
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that the vapour tension is equal to the partial water 
vapour tension of the air. The driving force has then 
become nil. 

Conclusions 
The humidity chart is an aid by which we may predict 
the course of the spray drying process, when only the 


100% 50% 5% 





| 
! 
3! 3:3 | 
; r +t 
ter t 
Fig. 12. Drying of a wet material (constant wet bulb 
temperature). 
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vapour tension curves of the material as functions of 
temperature and humidity have been measured. This 
has already been done for some materials, for instance 
for powdered egg by GREENE et al. [8]. 

The process can be accurately worked out by means 
of the vapour tensions, according to the example of 


H 




















% t 
Fig. 13. Efficiencies of 3 processes: ZL 1, ZL 2 and ZL 3. 


ty ty 


Fig. 11, if the temperature and humidity of the in- 
coming air are known. At each point of the process 
the driving force J H is then directly read off, and then, 
with the aid of the equation already derived, viz. 


dX 46.1077 1H 


oe yd 
dl ga \t + O23 Np, | 


(9) 
the fractional evaporation velocity per metre of height 
of fall of the drop at this point can be calculated. The 
diameter of the drop must be known and it must be 
remembered that (9) is the evaporated fraction of the 
whole drop when this has the same density as water. 

The above eq. (9) is correct, also for the evaporation 
of wet, drying spheres, if the rate of drying is con- 
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trolled by the resistance to diffusion in the air film 
round the drop. It is not until the moisture transfer 
within the drop starts to control the rate of evapora- 
tion that (9) loses its validity. 


Drying at constant wet bulb temperature 


Another conclusion is concerning the fact (depicted 
in Fig. 11) that during drying the wet bulb temperature 
rises according as drying proceeds. When a thermo- 
labile material is dried care will have to be taken that 
the highest wet bulb temperature (which the dry 
particles have at the end of the process) remains below 
a certain critical temperature. This means that at the 
start of the process the wet bulb temperature is un- 
avoidably low, so that the driving force is small and 
the process proceeds too slowly. It would be ideal if 
the wet bulb temperature could be kept constant just 
below the critical temperature. This can be achieved 
by starting the process with very hot air, with the 
wet bulb temperature at its maximum. Then this 
temperature is kept constant artificially at the highest 
possible level, i.e. is prevented from rising by con- 
tinuously admixing cold air along the entire length 
of the plant. The air point then does not intersect 
the adiabatic line, but shifts rapidly to the left, in 
such a manner that the wet bulb temperature remains 
constant (see Fig. 12). Apparently, therefore, parti- 
cularly at the beginning of the process, advantage can 
be taken of a higher driving force than by the known 
(adiabatic) process; cf. AH with AH’ in Fig. 12. 


8. THermat Errictency anp Erriciency Ratio 
We shall first deal with the evaporation of pure water 
and consider this process in the humidity chart, see 
Fig.13. Let us assume that the condition of the outer 
air is given by Z, and that after passing of the air heater 
condition L is reached. We shall consider three dif- 
ferent processes : 

(1) adiabatic process of infinite duration: ZL 1, 
(2) adiabatic process of finite duration: ZL 2, 
(3) non-adiabatic process of finite duration: ZL 3. 
In (J), the air is entirely saturated with water 
vapour, so that the final temperature ¢, is reached. 
This is the ideal process, with maximum evaporation. 
In (2), the air point passes along the adiabatic 
cooling line, without complete saturation being achiev- 


ed. Final temperature f,. 


In (3), heat is lost by poor insulation of the plant. 
Consequently the final temperature ft; is lower than ty, 
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though the evaporation performance is equal to that 
of process (2). 
The thermal efficiency is defined as 
heat of evaporation (28) 
heat input 
As the specific heat of the wet air is practically con- 
stant *, the ratio of the heats can be replaced by the 
ratio of the temperature changes, so that the following 
values are found in the three above-mentioned pro- 


(2) and (3) 


These thermal efficiencies are still dependent on 
the non-essential temperature of the outer air f,; there- 
fore, they are not a correct indication of the quality of 
the process. It is possible to obtain an expression 
(efficiency ratio) which does not contain this tempe- 
rature f, by referring the thermal efficiencies to the 
ideal process (/): 


(/) 9, = 1. (2) and (3) "r 


4 
Process (3) differs from (2) in that the effluent air is 
colder and in that it requires a larger apparatus. The 
latter fact is evident if attention is paid to the driving 
force AH’, which has become smaller than 4H owing 
to the lower wet bulb point P’, so that evaporation also 
proceeds more slowly. 

It is clear that efficiencies are particularly high 
if (other temperatures being constant) the temperature 
of the inflowing hot air t, is kept as high as possible. 
This possibility is limited by the increasing difficulty 
and expense of air heating and by the highest per- 
missible temperature of the drying drop. For, as the 
drop dries and the vapour tension becomes lower in 
consequence, the drop temperature will approach more 
and more closely that of the air passed off (ef. Fig. 11) 


Recirculation of spent air 


Recirculation of spent air makes it possible to increase 
the thermal efficiency by other means than those 
described above. This is shown in Fig. 14, where the 
same evaporation performance of the normal process 
Z 1, E is obtained by the recirculation process Z'L’'E’. 
The outer air Z is mixed with spent air EZ’, forming 

* To be exact, the spec. heat is (1 + 2H) kJ/kg °C for 
air with an absolute humidity of H. 


+ The useful evaporation heat is equal in processes (2) 


ad (3). 


: Spray drying 


the mixture Z’. When the two processes in Fig. 14 
are compared it is found that the thermal efficiency 
has actually risen: 
ty 
—t’ 


All air humidities have increased considerably, 
however, and the driving force AH has therefore 
decreased everywhere, so that the length of the ap- 


paratus is to be enlarged. 


Another disadvantage of recirculation of spent 
air is that the ventilation energy required soon be- 
comes too large, making the plant uneconomical. 


H 




















te 
Fig. 14. Efficiencies of normal process: ZLE, and of 
recirculation process: Z L’ E’. 


Drying of solutions and suspensions 


It was already shown in section 7 that when a suspens- 
ion or solution is dried the wet bulb temperature is 
not constant, but rises, approaching the air tem- 
perature in a process of infinite length (see Fig. 11). 


In an ideal process the position of this final tem- 
perature t, depends on the position of the vapour tension 
curves of the dried material and on the ratio of air and 
liquid in the apparatus: the more air, the higher is t,, 
and the lower the final moisture content of the mate- 
rial. 

If we compare this case with that for pure water 
(see Fig. 10) the situation is found to be distinctly 
less favourable: at equal H levels the driving force 
AH is greatly reduced; it is even impossible to get 
beyond H.. It is clear that the evaporation perfor- 
mance of the same plant has therefore become con- 
siderably lower in comparison with the evaporation 
of pure water. If these performances have to be 
compared, however, it is necessary to take into con- 
sideration the difference in driving force 4H. This 
can be done with the aid of the transfer unit concept 
introduced by Catton and CoLBurn [9], as explained 
in the next section. 
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9. TRaNnsFEeR Units 
These dimensionless numbers were introduced by 
CuILTon and CoLaurn [9] for processes of heat transfer 
and mass transfer (distillation, absorption, heat ex- 
change, etc.). The difficulty of a diffusion process may 
be characterized by the number of transfer units, 4). 
For instance, in heat transfer through an air film .N, is 


defined by: 


3 “dt - 
N, | it (=) 


Let us first consider evaporating drops of water 
in air of which the surface temperature is the wet bulb 
temperature - adiabatic saturation temperature. 
With adiabatic air wetting (¢/. Fig. 13) from L to 2 


a 


N; 


‘ 
“dt t,—t, 
| It t,—t 


é 


(30) 


This simple relation [eq. (30)} between efficiency 
ratio and the number of transfer units only applies for 
adiabatic evaporation of a pure liquid. 

Next, let us consider the case of the drying of 
drops of a solution or suspension. We have already 
seen (Fig. 11) that the vapour tension curves are lower 
in the humidity chart as the humidity is lower. It is 
now possible to draw full advantage from the special 
The defin- 


ition eq. (29) remains perfectly valid, while At can 


properties of the transfer unit concept. 


be read off in the humidity chart step by step. As 
At decreases more rapidly, the number of transfer 
units will be much larger than in an equivalent evap- 
oration (i.e. with equal increase of H) of pure water. 
A larger apparatus is therefore also needed. 

Apparently the drying of a particular solution in 
a spray drier should be evaluated as follows 

1. Using pure water, the thermal efficiencies and 
the number of transfer units are determined as de- 
scribed in section 8 and with the aid of the above 
eq. (30). 

2. The vapour tension curves of the material to 
be dried are determined as functions of temperature 
and humidity. If necessary, the maximum tempera- 
tures permissible to avoid thermal decomposition of the 
material are also determined. 

3. The course of the process is plotted in the humid. 
ity chart. The ratio of air and liquid currents is 
adjusted so as to obtain the desired moist content of 


Spray drying Engincering Science 
the product required with the same number of transfer 
units as determined under 1. Attention is also paid 
to the permissible temperature of dry air passed in. 

When the drop distribution during the spraying of 
the solution is the same as during the spraying of pure 
water, the measured number of transfer units of the 
apparatus remains the same. When this is not the 
case, allowance must be made for the influence of 


particle size (see below). 


Example of calculation 
Water is sprayed in a spray drier. The following temperatures 
are measured: 
a) outer air is 15° C (59° F) dry bulb thermometer, 
b) outer air is 7° C (44° F) wet bulb thermometer, 
¢) air after heater 115)° C (240° F) dry bulb thermometer, 
d) discharged air 57° C (135° F) dry bulb thermometer, 
e) discharged air 32° C (89° F) wet bulb thermometer. 
The process can be easily followed on the humidity chart 
Fig. 9. 
115) 73) 
115) 15 
temperature of 73)° C (164° F) having been found as follows. 


1. The thermal efficiency is 9 42%, the 


From a) and b) it foliows that the absolute humidity is H = 
0-003. 0-003 and t= 115)°C 
[240° F]) one follows a certain adiabatic cooling line, which 
ends at the wet bulb temperature of 35° C (95° F). But the 


From the air point (H 


measured wet bulb temperature of the spent air is 32° C 
(891° F). In combination with the dry bulb temperature of 
57° C (135° F) this points to an absolute humidity of H 
0-020 by the 
adiabatic cooling line and find the temperature required to be 
73.° C (164° PF). 


0-020. We intersect the horizontal line 


»o 
- 0° 


80% 
value of the thermal efficiency of 80% having been found as 


2. The efficiency ratio is n, 525%, the maximum 


follows. The above-mentioned adiabatic cooling line pointed 
to the wet bulb point of 35° C (95° F). The required maximum 
value is therefore 

115) — 35 


: - == 60%. 
115) 15 


3. The heat loss follows from the fact that in an adiabatic 
process the final temperature, according to 1., might have 
been 731° C (164° F). In fact, the temperature of the spent 
F), so that 73) — 57 
Calculated on the total drop in temperature from 115)° C 
73), — 57 
115 — 57 


4. The number of transfer units. 


air is 57°C (135 16)° C was lost. 


this is 284%. 
Here we have a non- 
adiabatic process, so eq. (30) may not be applied. We therefore 
use the definition eq. (29): 


115) — 734 


log mean Jt ° 


(31) 


n= {4 
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In (31) only the useful (evaporation) heat is, of course, consid- 
ered, so that the drop in temperature to 73° C is taken (not 


to 57°C). The driving force varies from (115) — 35) to 
(57-32). The logarithm mean of this is 46° C, 


ql 79) 
Therefore, N, 115} 734 


é 0-91 transfer unit. 


Transfer units on a basis of concentration 


It was said at the beginning of this section that the 
transfer units can also be calculated on a mass transfer 
basis. For spray drying the following definition is 
obvious: 
’ “dH ee 
“tm | 1H (32) 
These transfer units are exactly proportional in 
number to those calculated on a temperature basis 
according to eq. (29) as 
Nim . (33) 
Proof 


wet bulb line: 


Use eq. (25) and (26) of section 6: 


h.At=k, iAH 


adiabatic cooling line: 
Ad. 
Hence, 
1H h 
It kim 4 


Theref« re 
[ dH 
1H h 


km Cp [ dt 


it’ 


from which follows the above-mentioned result, 
eq. (33). 
For the special case in which Lewts’ equation 


applies, see eq. (27), k,, C, — h, so that then 


N, N, (36) 


“"tm 


This obtains for the combination water/air, and for 
drying aqueous suspensions or solutions, regardless 
whether the adiabatic cooling line is straight or 
curved. 


Calculation of the height of a transfer unit for falling 
drops 

We shall now calculate the height required for a pro- 
cess of one transfer unit, for water drops falling in 
stagnant air, starting from the equation for fractional 
evaporation of a falling water drop per metre height of 


1H 


= 1 + 0-23 Np’). 


46.107 (37) 
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Suppose the flow of air through the tower to be 
V kg/see and the flow of liquid at a given point 
Lkg/sec. At that point the humidity of the air rises by 


70-5 


dH L ditw 4 
, tw _ 4.6.10-7 © 44 (1 1 0.93 N2*). (38) 


dl , V oP 


Now, the height equivalent to one transfer unit is 
defined as 

length dl 
; (39) 
number of transfer units in dl 


2.2 -10° +0* 
=) 
*, 


(ob: 
O23 N 


en 


a“ 





o2 03 O4 ° °) o8 OF ae 


Fig. 15. Height equivalent to a transfer unit. Free falling 
water drops in stationary air (0 1 kg/m*). 


From (38) we can find by means of (39) the value of 
H,, as follows: 


dl dl 22.18°0R V 
l 


4) 
(Mjar dH/AH +oas nes: ™ (4) 


H, 
where L = kg water/sec at the given point. 

Eq. (40) therefore obtains for every moving drop 
of water in air; hence also for drops falling at a con- 
stant settling speed. This speed was derived in sec- 
tion 3. When these speeds are substituted for water, 
the value of H, for drop diameters between 0-1 and 
1-0 mm follows as shown in Fig. 15. 


This result can be interpreted as follows. If the 
job of humidifying air by a water spray is considered, 
it appears that the difficulty of this job (H,) is in 
direct proportion to the ratio of air and liquid (V/L). 
Further, the smaller the drop size the more this dif- 
ficulty decreases (Fig. 15). 

Next we will treat H, for falling drops of a solution 
or suspension. We suppose that the feed J at a given 
moment only contains a fraction B of water and 1-B 
of solid matter. For the fractional evaporation of the 
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falling drop we then get the modified equation 


46.107° AH 


dx, {] 
dl B “oP 


« 0S 
23 J He } 


(37°) 
because the evaporation is taken to be equal to that 
of a drop of pure water. 

At that moment the humidity of the air rises by 


i at BL... ee 
, 0-7 
dl yo dl y 46.1 » De 


vO 5 
(1 + 0-23 Np.”) | 


(38) 


which therefore tallies entirely with eq. (38) if L is 
taken to be kg feed/sec. The same also applies to H,, 
represented by: 


2-2. 10%» 1p 
H, 1 + 0-23 NS? ° 


(40) 


where again L ku feed sec at the given point. 

The result, viz. that the H, is the same regardless 
of the solids content of the drop, may be interpreted 
as follows: for the wetting of air it makes no difference 
at all whether the drop contains less or more solid 
matter (a difference in vapour tension, 4//, is auto- 
matically accounted for in the concept of transfer 
units). On the other hand, for the drying of the drop 
it is of course advantageous if the solids content of the 
drop fed in is as high as possible. 

Analogous considerations obtain in the drying of 
hollow drops as compared with solid drops of the same 


size (see section 5). 


10. HyproGen anp Steam as Dryine Gases 
Although air is an obvious drying gas, there may be 
some other gas that has more advantages for this 
purpose. 

Hydrogen possesses the advantage of having a 
higher rate of diffusion and heat conductivity. Care- 
ful analysis shows that the rate of evaporation of 
water in hydrogen is indeed three to six times as great 
as that in air (depending on the case under considera- 
tion). There are, however, considerable practical 
disadvantages, such as expense and inflammability, 
that make its use impracticable. 

Superheated steam can also be employed as a drying 
gas. It is interesting to note that the evaporation 
rate of wet material in a steam atmosphere is only 
determined by the heat transfer of the drops. In the 
case of evaporation in air the material diffusion 
coefficient is also a determining factor. In the steam 
process there is, of course, no question of material 
diffusion in the gas phase. 


Spray drying 


Re hydrogen 
The influence of the drying properties on the rate of 
evaporation is evident from FrOssLiInGs evaporation 


equation: 


=~ x05 yy0-33) kmol 
(1 + 0-276 Np, Ng, ) —_* 

The rate of evaporation depends on the following 

factors: 


1. diffusion coefficient, 

2. difference in concentration, 

3. ReyNoLps’ number and 

4. Scumipts number. 

1. The diffusion coefficient D,. This depends equally on 
the magnitude of the molecular weight of the evaporating 
material and of the gas. This is apparent from the equation 


(on an empirical basis) for molecules A and B: 


7 
» LP | 
+ (Vin) B 


M, | Mp (41) 


P (Van) 4 


Further particulars on this correlation eq. (41) are given by 
GILLILAND [10}. 

If we compare the rate of diffusion of water (M = 18, 

Vy, = 1S) in air (M = 29, V,, = 29-9) with that of water in 

14-33), we find the value of D, in 


hydrogen to be 3-2 times as high as in air. 


hydrogen (M = 2, V,, 
No other gas has such a favourable rate of diffusion. 
10, (kmol/m*) 


on either side of the gasfilm. This quantity depends on the 


2. The difference in concentration A p, RT 


difference in wet bulb and dry bulb temperature, as discussed 
in section 6. At a given gas temperature we now compare the 
wet bulb temperatures in hydrogen and in air; ¢, is found 
from the gas point in the humidity chart by drawing a line 


with a slope according to eq. (25) 


1H h 


‘ _ } (42) 


Now we can derive &,, and A from Fréss~ines equation for 


mass transfer [eq.(2)] and its analogon for heat transfer 
feq. (4)]: 
? (I 


+ 0276 NPS NE*), 


> (43) 


2 kp l 0-276 vos you 
Dp Ut OS Nhe “Pr )- 


(44) 
For low values of Np,, when 1 > 0-276 Np NS, we find 


accordingly : 


1H 


PP (45) 


and for high values of Np,, when 1 << 0-276 Np? NS: 


1H h Up Noo. 
At ky 3 , “te 


(46) 
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With eq. (45) and (46)* we may estimate how the slope 1H/11 
will be affected, when hydrogen is used instead of air. For 
low values of N,,, it follows from eq. (45) and the physical 
properties that 


Awater 


(CO, Nadu (145.21) (6). Ny jai: 
Therefore, the slope 14/11 will be 30 times steeper than in 
the case of air. For high values of Np, , it follows from eq. (46) 
and the physical data that 


ve 67 


7 ae 


l “ 
— (14-5. (2-1)°°) (0 
4. 


Dl a 
3 (f ” Nie), 


Therefore, the slope will be 24 times steeper than for air. 

It would now be possible to construct a new humidity 
chart for the combination water/ hydrogen, for instance by 
the method indicated by Barta and Garper [11]. In this 
particular case, however, one can still make use of the known 
chart for water/air. For the humidity seale must be magni- 
fied 14-5 
greater than that of H,. At the same time the slope of the 


<, in view of the fact that the density of air is 14-5 » 


adiabatic cooling line [eq. (26)] in this chart: 


du c 


dt } (47) 


14-5 (€C 


Consequently, one can for practical purposes use the air 


is also 14-5 = greater, as (C’,) i. pair’ 

water chart for the case of H,/water, including the adiabatic 
lines. The above-derived slopes for the wet bulb point, 4H/At 
should also be divided by the same factor 14-5 before use in 
the air/water chart. We then find, that the slope should be 
multiplied by 2-1 for low values of Np, and by 1-64 for high 
values of N,.,. 

Example. Suppose N,, has a low value. Suppose the 
gas point is HW = 0, t= 200° F. For 
wet bulb point: H = 0-026, t = 85°F. Therefore the driving 
force JH = 0-026. For hydrogen, we 
0-026 
200-85 
103° F. So the driving force is 1-8 times 


air we find in Fig. 8: 


construct a wet bulb 


point line with a slope of 2-1. . The wet bulb point 


is H = 0-047, t 
larger than with air. Similarly, for a high value of V,, we 
find a driving force 1-46 times larger than with air. 

3. Reynoips’ number Np, = vD oju. When we consider 
falling drops, the constant settling velocity will, of course, be 
much higher in hydrogen than in air. Let us proceed with the 
case that Stokes’ law obtains (that is for low N,,). Then the 


* If we try to check these equations for the case of water 
, } , , ’ 
and air, we find for low Np,: i a N,, and for high N»p,: 
h Rai 4 , 
As was mentioned in section 6, eq. (27), it is 
kn Cp 


an experimental fact that 


- Nz". 
~ 1 for water and air. On 
in Cp 

the other hand, N,, = 0-917 according to Jakons “Heat 


Transfer,” so that Ni’ == 0-04. 
psychrometric ratio” 


The equation for the 


rT — nui 
m\p 
checks quite well with a graph in GriuiLanp [12]. 
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practically unchanged drop weight is 


, SC , 
3au De 4 IP gq Ao. (48) 


Therefore » - 
nos — (29 Ao | 


Ke Isp? (49) 


So, comparing hydrogen with air, in Stokes’ case Nee is 
0-069 oats ‘ : 
— 0-52 times larger. A comparable factor will hold in 

“52 


the case of high values of Np, ' 


| When the 
o Dp 


substituted for water in hydrogen we find Ng, 


4. Scumiprs number Ny values are 


1-43, as 
compared with 0-60 for water in air. 

Summing up the four factors just discussed, we 
find that, according to Frésstings eq. (2), when 
changing over from air to hydrogen as drying gas the 
following changes occur: 

(a) At low values of Np, 
according to the first term of eq. (2), in such a way 
that the rate is 3-2.1-8 

(b) At high values of Np,: according to the second 


: evaporation is effected 
5-8 times higher. 


term of the FréssLivec equation, the rate is increased 
2 0-33 

| = 928. 

0-60 | 

The required height of fall is affected less favour- 


by a factor 3-2. 1-46.0-52. 


ably, because the constant settling velocity in hydro- 
gen will be higher than in air, as already mentioned 
above. 

When hydrogen is used instead of air as a drying 
medium, a number of technical difficulties may be 


Table 1 





Air H ydrogen Steam 


Dyn. Viscosity at SOC, mw, 
1-2. 1075 


18 


12.107° 
14-4 


js eee ee 
C,, in kd/kg © 
conductivity, &, in 


Spec. heat, 
Heat 
J/m. see °C » is 0-03 0-21 0-02 
Molecular weight, M, in 
ND a 6a aS ae ee 29 - 18 
Diffusion coefficient with 
H,0, D, in m*/see 


Molecular volume JV, [10] 


3.10% 9).10°5* 
29-9 14-3 15 

Thermal diffusivity in 

18.1075 12-9.10°° 5-2.10°° 
0-74 0-74 0-78 


m?/sec = 
Number of Pranpri 
Number of Scumrpt 
(with H,O). 
Number of Lewis 
(with H,O). 


0-60 


0-917 





above. 


* By calculation, see 
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expected, as, for instance, the recirculation of this 
comparatively expensive gas, necessitating the appa- 
Although the 


use of hydrogen as a drying gas has been suggested by 


ratus to be gas-tight and fire-proof. 


some authors [12], [13], so far these technical dif- 
ficulties have prevented the realization of H,-drying. 
Some physical constants are cited of air, hydrogen 


and steam. 


Re superheated steam 
Another drying gas previously suggested [15], [16] 


is superheated steam. When a closed recirculation 


PRESSURE WATER WET MATERIAL 


PRESSURE OF STEAM 








= 





aa | 
=a 


Fig. 16. Wet bulb temperature in superheated steam. 


system is used, the consumption of heat will be a min- 
imum, as only the evaporation heat of water has to 
be supplied (apart from heat losses). 

Of course, in evaporation in overheated steam 
there is in the gas phase no material diffusion, but only 
heat transfer as a limiting factor. The surface tem- 
perature of a drop is determined entirely by the steam 
pressure (independent of the steam temperature): 
1 ata corresponds to 100° C, } ata to 81°C, ete. (to 
be found in the humidity chart or a vapour tension 
table). 
lower equilibrium vapour tension. Conversely, there- 


At equal temperatures wet material has a 


fore, at a given steam pressure, the wet bulb tem- 
perature of wet material will be higher than of pure 
water. This is illustrated in Fig. 16. 

If one wishes to compare the rate of evaporation 
with that in air, this can only be done if a reasonable 
basis of comparison is agreed on. This may be done 
as follows. 


Drying in air—Using dry air of temperature ¢, the 
corresponding wet bulb temperature is ¢,,; the driving 
force is then t, —ty. 


Drying in superheated steam—Such a pressure is chosen 
that the corresponding temperature (according to the 
vapour pressure table of saturated steam) is ¢,,.. Further 
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the temperature of the steam is chosen ¢,; again the 
driving forces is t, —t,. 

In order to compare the rate of heat transfer, 
eq. (44) should be used. 


of Np,, only the first term is of interest. The value of 


In the case of low values 


the heat conductivity, k, is for steam 2/3 times that 
for air, so that the rate of heat transfer is also de- 
creased by a factor 2 3. 

In the case of high values of Np, , only the second 
Now Np, and Np, also affect h. 
For steam, Np,’ is nearly equal to the value for air. 
But Np, should be calculated for a constant settling 


term is of interest. 


velocity in steam. Then we find 


2 | Np, in steam }05 


hee hai ; ; 50 
neoauen =" 3 Nye in air (50) 

As in the case of hydrogen, the practical difficulties 
entailed by the recirculation of the gas are great. 
Moreover, the removal of the product from the steam 
atmosphere, without condensation of water, is an- 


other practical difficulty to be faced (see [15}). 
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NOTATION 
Surface area, m? 
Fraction by weight of water in feed 
Drag coefficient for sphere in fluid, dimensionless 
Specific heat of gas at constant pressure, Joule/ 
kg °C 
Vapour concentration in air, kmol/m* 
Diffusivity, m*/sec 
Diameter of droplet, m 
Energy, Joule 
Extra evaporation in eq. (22) 
Force, Newton 
Acceleration of gravity, 9-8 m/sec* 
Absolute humidity, kg water vapour/kg dry air 
Height of a transfer unit, m 
Absolute humidity for a given wet bulb tempera- 
ture 
Final humidity, reached after infinite time 
Coefficient of heat transfer, Joule/m*sec °C 
Heat conductivity Joule/msec °C 
Coefficient of mass transfer, kg/m*sec (unit of 4H) 
Length, height or path, m 
= Liquid rate, kg/sec 
= Molecular weight, kg/mol 
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Mass, kg 
Rate of transfer, kmoles/sec 
Number of NusseEt, dimensionless 
Number of Pranpti, dimensionles 
Number of Lewts, dimensionless 
Number of Reyno.ps, dimensionless 
Number of Scumrpt, dimensionless 
Number of transfer units, dimensionless 
Number of transfer units (mass basis), dimension- 
less 
Partial pressure, N/m? 
Total pressure, V/m? 
Flow of heat, Joule/sec 
Gas constant, 8-3 . 10° Joule/kmol °K 
Absolute temperature, ° K 
Temperature, °C 
Temperature of air, °C 
Wet bulb temperature, °C 
Final temperature reached after infinite time, °C 
Vapour rate, kg/sec 
Velocity, m/sec 
- Molecular volume acc. Korr-Le Bas, see eq. (41) 
Weight of liquid, kg 
- Fraction by weight, dimensionless 
Finite difference 
Thickness of film, m 
= Thermal efficiency, dimensionless 
: Efficiency ratio, defined in section 8. 
Dimensionless time, defined in eq. (16) 
Heat of evaporation, Joule/kg 
su = Dynamic viscosity, kg/m sec 


- Kinematic viscosity, m?/sec 
= Density, kg/m* 
Density of liquid phase, kg/m* 
Density of gas phase, kg/m* 
Density of sphere, kg/m* 
Dimensionless path, defined in eq. (17) 


- Time, sec 
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The catalytic polymerisation of linseed oil with sulphur dioxide 
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(Received 14 November 1951) 


Summary —The catalytic polymerisation of linseed oil with sulphur dioxide was investigated. The presence 
of small amounts of oxygen proved to be indispensable for this reaction, complete absence of oxygen resulting 
in an activation (conjugation) of oil without a marked increase in viscosity. Stand oil was continuously 
prepared on a small scale in a column at a relatively low temperature (approx. 300° C), using sulphur 


dioxide as a catalyst. 


Résumé 


On a examiné la polymérisation catalytique d’huile de lin avec du bioxyde de soufre. La présence 


de faibles quantités d’oxygéne parut nécessaire A cette réaction parce qu'une absence totale de l’oxygéne 


résulte & une activation (conjugaison) d’huile sans une augmentation considérable de la viscosité. 
La préparation continue d’huile siccative fut effectuée a petite échelle dans un appareil a colonnes a 
une température relativement basse (300° C environ) en utilisant du bioxyde de soufre comme catalyseur. 


INTRODUCTION 
The polymerisation of linseed oil with sulphur dioxide 
as a catalyst has been the subject of extensive studies 
in this laboratory. The catalytic action of sulphur di- 
oxide was first described by WATERMAN and van VLo- 
pRoP [1] and was also applied on a technical scale. In 
later experiments [2] preparation of stand oil was per- 


10 Chem, Eng. Sci. Vol. 1 


formed continuously, using relatively high tempera- 
tures (360° C). 


* This research work was carried out by the collaboration 
of the Chemical Engineering group (Prof. H. I. WaTeRMAN) 
and the Mechanical Engineering group (Prof. E. F. Boon) 
of the Department of Chemical Engineering, The University, 
Delft, Holland. 
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The action of the catalyst is based on an activation 
of the linseed oil molecules, bringing about a con- 
jugation of the linoleic and linolenic acid groups, 
which, according to KaprrELMEIER [3], polymerise 
easily at higher temperatures, giving cyclic reaction 
products [4]. 

In accordance with this conception, the polymeri- 


sation of linseed oil activated beforehand (e.g. with 
sulphur dioxide at 180°-200°C and 60 atm [5] or 
with other catalysts [6] takes place rather easily; 


irons ~ 
Vibrator 


the polymerisation of tung 
oil, consisting for the greater 


part of conjugated acid gly- 





cerides, is not greatly af- 


fected by sulphur dioxide. 
50 


r In the catalytic manu- 
ong Air 


facture of stand oil on a 
larger scale, difficulties arose 
in the mixing of oil and 


catalyst. A bell-stirrer, satis- 





factory in laboratory experi- 
to bad 
in a semi-technical plant. 


ments, led results 


Therefore another system 








was used for continuous 








stand oil preparation, using 
a packed column as_ the 
Fig. 1. Batch polymerisa- reactor. 


tion of linseed oil. As the influence of oxy- 


gen (air) on the polyme- 
risation process cannot be neglected (this, of course, 
may be important when working in a closed appara- 
tus), this problem was studied more closely in the 


following experiments. 


Batcu EXPERIMENTS 

Batch polymerisation experiments were made in a 
Pyrex reaction tube (contents approx. 1} 1), using 
400 ml of a commercial bleached linseed oil in each 
experiment (Fig. 1). The tube was shut by means of 
a cork and heated in a metal bath. Stirring was 
performed by means of a vibrator [7], working as 
a churn and causing an intimate contact between oil 
and gas. Contact with the air was prevented as much 
as possible by using a fairly powerful stream of the 
gaseous catalyst and a narrow outlet-opening. 

Two series of experiments were carried out, using 
different mixtures of sulphur dioxide and oxygen, 
at 300°C (45 min) and 330°C (15 min) respectively. 
The oil was heated to the desired temperature in 
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10-15 min, using carbon dioxide as an inert gas, 
which will have caused a fairly complete removal of 
peroxides present [5]; after the polymerisation the 
liquid was similarly cooled in a stream of carbon 
dioxide. The refractive indices (nj), the viscosities (7) 
and the acid values of the reaction products were 
determined; the results are collected in Table 1 and 
Fig. 2. 


Table 1. Batch polymerisation of linseed oil with sulphur dioxide 
and air. Amount of oil 400 ml in each experiment; oil preheated 


to the desired temperature in a CO,-atmosphere 





A) Temp. 300° C, 


reaction time 45 min 


B) Temp. 330° C, 


Catal yst : , # : 
y reaction time 15 min 


1) acid 20 1) acid 


nD (poise) value "D (poise) value 
1-4902 Is 
1-4902 
1-4910 
1-4909 
1-4896 


1-4854 


1-4885 16 
1-4891 18 
14905 21 
1-4905 32 
1-4894 18 
1-4837 15 
Linseed oil 14814 


(original) 





75 min 330°C 
% min 300°C 


50 
Air in catalyst 
Fig. 2. Batch polymerisation of linseed oil 
with SO, and air. 


It will be seen that air does influence the catalytic 
action of the sulphur dioxide, the most favourable 
the 
30 vol.% of air. The total amount of oxygen used 


composition of catalyst containing approx. 
in the experiments was comparatively low; oxidation 
of the oil was therefore of little importance. 

In some other experiments efficient precautions 
were taken to avoid any influence of traces of oxygen 
in the polymerisation experiments. The linseed oil 
was pre-treated with nitrogen (15min at 290° C) 
in order to destroy any peroxides present in the 
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‘able 2. Batch treatment of linseed oil with sulphur dioxide in the absence of oxygen. Amount of oil 50 ml. 





Properties of reaction products 


Diene No. 


Reaction conditions E 
“LLIS 


Viscosity 


poises 


I. Vv. Wiss a id val uc 
and 


JONES 


0-3063 
0-3066 
0-3058 
03060 


1-4808 
1-4861 
1-4879 
1-4874 


A. Original oil 
l hr 285° C, 
lhr 290°C, 
lhr 200°C, 


0-9288 
0-0365 
0-9419 
0-9405 


8Vhr SO, 
8 hr } 
12 Vhr 


0-9283 
0-9363 
00-9404 
0-939] 
0-9474 


03065 
0-3062 
0-3057 
0-3063 
3047 


1-4808 
14852 
1-4867 
1-4871 
1-489] 


B.* Original oil. . . . 
Lhr 300° C, 7 Vhr 
| hr 300° C, 7 Vhr 
‘hr 300°C, 7 Vhr 
lhr 300°C, 7 Wher 





*In these experiments the sulphur dioxide was not treated with pyrogallol. 


oil [5]; also, in some cases, the sulphur dioxide was 
distilled from solid pyrogallol [8]. 

These precautions had a remarkable influence on 
the results of the experiments, as is shown in Table 2. 
Although the refractive indices of the reaction pro- 
ducts showed a considerable increase, the viscosity 
remained very low, and also the specific refrac- | 
tion did not alter very much. The diene values, on 
the other hand, were rather high. These observ- 
ations may be compared with the results of 
the experiments of WATERMAN c.s. [5], in which 


linseed oil was “activated” by heating with 
sulphur dioxide at a temperature of about 
200°C under a high SO, pressure, the activation 
being more marked, however, in the latter ex- 
periments. 

From these results it may be concluded, 
in accordance with the views of WATERMAN, 
c. 8. [1], [9] that the catalytic influence of 
the sulphur dioxide is restricted to the con 
jugation of the linseed oil fatty acid groups. 
Activation of linseed oil appears possible not 
only with high pressure sulphur dioxide, but also at 
atmospheric pressure, albeit at higher temperature. 
For the polymerisation of linseed oil small amounts 
of oxygen are indispensable. It may be advisable, 
therefore, to add a small percentage of oxygen to the 
catalyst in the catalytic (SO,) manufacture of stand 
oil; in most cases, however, the oil to be heat bodied 
will contain a certain amount of oxygen as peroxides, 
which is sufficient for a normal course of the poly- 
merisation process (compare Table 3). 


10* 


CONTINUOUS POLYMERISATION OF LINSEED OL 


Continuous polymerisation of linseed oil with sulphur 
dioxide was performed in a packed column (40” » 1’) 
of stainless steel, filled with aluminium turnings and 
externally heated by means of a Dowtherm system. 
Oil and gas were introduced into the bottom of the 





| - t Vacuum 
Packed G 
. 


Column | 














— 





=) Dowtherm system 
Fig. 3. Continuous polymerisation of linseed oil. 
column and passed the reactor in concurrent flow 
(see Fig. 3). The oil was pre-heated in a metal bath 
D, the oil flow could be measured in B and was 
metered by means of a small gear pump C provided 
with a bypass valve R. The sulphur dioxide and air 
flow could be varied by means of the capillary tubes 
K and L and the heights of the mercury columns in 
the pressure regulators M; the gases were conveyed 
via the flow meters J in a mixing vessel H and intro- 


The oil and 


duced into the bottom of the column. 
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Table 3. Continuous polymerisation of linseed oil (non pretreated ) with sulphur dioxide 


and air in a packed column. Volume of oil in the column 250 ml 
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oil apparently being sufficient for 
the polymerisation (expts. 2 and 5, 





Catal yst 


Expt. Temp. Linseed Reaction m 
No, c SO, air ou oom “D 
lar lhe ml min mon 
I 200) l 0 8 43 
2 200~— I 0 0 OO 1-4900 
3 200 l 0 30 64 1-4905 
4 200 l 0 340) 83 14910 
5 200 7 O-3 4-05 rh) 1-4889 
6 200 7 O-3 tS 66 14809 
7 200 «7 O3 31 sl 1-4903 
8 300 (OT O-3 50 OO 1-48908 
of) 3100 OF 0-3 765 33 14900 
10 32000 «7 O-3 Sl 31 1-4913 
ll 3000 «O75 0 48 52 1-4907 
12 300 10 0 5 Mw 1-4907 
13 300 1-25 0 5 ww 14910 
14 3000-50 0 49 51 1-4910 
- Linseed oil 1-4808 


Viscosity 


3 and 6, 4 and 7, 8 and 12). 


in poises Acid The accumulation of free fatty 
Fa sl ee acids is one of the main problems 
in continuous heat bodying of oils. 
Because of the relatively low tem- 
13 I4 . 2 : 
198 “ peratures which can be applied, the 
26 14 acid values of the polymerised oils, 
33 16° although rather high, can be kept 
16 4 within the Dutch standard specific- 
1) 15 ations for zine white stand oils {11}. 
28 17 When using temperatures above 
26 18 300°C the formation of free fatty 
22 20 acids greatly increases, as follows 
33 24 from the experiments 8-10. 
- - The results of the continuous 
“ heat bodying process of linseed oil, 
99 17 as regards the reaction time for pre- 
paration of stand oils of a certain 
0-5 4 


viscosity, are in good agreement 





gas flow were adjusted in such a manner that good 
dispersion of gas in liquid was obtained; as was 
proved by preliminary experiments in a glass column 
small gas bubbles moved slowly upward in a regular 
stream, ensuring intimate contact between oil and 
catalyst. The reaction products, leaving the top of 
the column as a foam, flowed via the degasser G in 
the collecting vessel Z. 

In some experiments the column was drained to 
measure the amount of oil in the reactor, in order to 
determine the contact time of the reaction and the 
rate of polymerisation. This amount appeared to be 
all but independent of the reaction conditions as used 
in our experiments, and was fixed at 250 ml of oil, 

The application of the Dowtherm heating system 
made it possible during the experiments to maintain 
the temperature (thermometers 7,, 7,, 7;) along the 
column constant within 2°C. 

Some of the results are collected in Table 3. 

They show that continuous preparation of stand 
oil is possible at relatively low temperatures, using 
sulphur dioxide as a catalyst. 

Varying the amount of catalyst (0-75-1-50 hr of 
SO,, expts. 11-14) had no influence on the rate of 
viscosity increase. The addition of air appeared to 
be of no importance in these continuous experiments, 
the amount of oxygen in the non pretreated linseed 


with batch experiments made at the 


same temperature and with the most favourable 


catalyst composition (Tables 1 and 3). 

The packed column can easily be developed on a 
larger scale and be used for heat bodying of all types 
of drying oils, as temperature, flow of oil, amount of 
catalyst and eventually pressure can be widely varied. 

The authors wish to express their gratitude to 
Ir. N. J. Stenstra, Mr. W.A. Kiaassen and Mr. 
J. A. BrexkMan, for doing part of the experimental 
work. 
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Summary —A synopsis is given of a recent theory of distillation of multicomponent mixtures in which the 
product of absorption factors is calculated from a summation equation. A constant ratio of volatilities is 
assumed. The calculation of a distillation at variable reflux can be reduced to one at constant molar reflux. 
The theory is applied to a distillation at constant reflux. The summation equation is solved either by 
successive calculation of the absorption factor product or algebraically or by graphical methods. The 
successive calculation is by far the quickest if only one reflux ratio is considered. The algebraic method 
involves the solution of two algebraic equations of higher order than two. It is shown that the multi- 
component mixture may be approximated by two binary mixtures and if this is done quadratic or even 
linear equations have to be solved only. Even so, however, the algebraic solution takes up much time. 
The graphical method is also based upon the approximation by binary mixtures. It is preferred to the 
algebraic method on account of its quickness and it is easily surveyable. Graphs are given from which all 
possible columns and reflux ratio’s which lead to the desired distillate and residue, can be found. 


Zusammenfassung 
gegeben. 


Ein AbriB einer neuen Theorie der Destillation von Mehrkomponentensystemen wird 
Das Produkt der Absorptionsfaktoren wird aus einer Summengleichung berechnet. Konstante 
Destillation mit Ricklauf wird auf die 
Destillation mit konstantem Ricklauf zurtickgefiihrt. Eine Destillation mit konstantem Riicklauf wird als 
Beispiel behandelt. Die Summengleichung wird gelést mittels sukzessiver Berechnung des Produktes der 
Absorptionsfaktoren. Eine algebraische Losung und graphische Methoden kommen gleichfalls zur Ver- 
Von diesen fiihrt die sukzessive Berechnungsweise weitaus am schnellsten zum Ziele falls nur 
ein Ricklaufverhaltnis bericksichtigt wird. Die algebraische Methode beruht auf der Lésung zweier Glei- 


relative Flichtigkeit wird angenommen. Die veranderlichem 


wendung. 


chungen hdheren Grades. Das Mehrkomponentensystem kann aber durch zwei Zweikomponentensysteme 
approximiert werden und das eréffnet die Méglichkeit, die algebraische Lésung auf die Lésung linearer 
Gleichungen zurickzufthren. Dennoch sind die Rechnungen verhaltnismaBig zeitraubend. Die graphische 
Methode beruht gleichfalls auf der Naherung durch Zweikomponentensysteme. Es wird dieser Methode 
wegen ihrer Ubersichtlichkeit und Kirze der Vorzug vor der algebraischen gegeben. Alle mdglichen 
Kolonnen und Riicklaufverhaltnisse mit welchen man das gewiinschte Destillat und Bodenprodukt bereiten 


kann, sind mittels einiger graphischer Darstellungen zu finden. 


INTRODUCTION 


During the past five years several papers have appeared 


on the subject of distillation of mixtures consisting of 
many components. Alternative analytical calculation 
Harperr [1], 
UnpERWoop [4], Murpock [2] and Smrruvysen [3]. 
They are applicable to a distillation in which the molar 


methods have been developed by 


reflux is constant and the volatilities of the compo- 


nents have a constant ratio. In these methods an 
algebraic equation of the Nth degree is arrived at, 
when the system contains N components and the 
calculation of the distillation is in principle reduced 
to the solution of two such algebraic equations, one 
for the stripping section and one for the top- 
section. 

An extremely elegant method for calculating mini- 
mum reflex ratio has been given by UNDERWooD but 
the calculation by these methods of a distillation at 
a finite reflux ratio, which is higher than the minimum 
reflux, has still remained a laborious procedure, in- 
volving a large number of numerical calculations. 

One of the present authors has developed a theory 
in which the absorption factors of an arbitrary refer- 


ence component are calculated on the different trays 
in a distillation column. The equation of conservation 
of heat is used to calculate these absorption factors 
and therefore, a variable reflux can be taken into 
account [6]. The theory is applicable to the most 
general type of distillation in which the relative 
volatilities may vary, heat loss may occur or even 
the 


components. In practical calculations, however, one 


chemical reactions may take place between 
has always first to consider the case of a constant 
relative volatility and further refinements are taken 
into account by successive approximation. No trial 
and error is involved thereby. 

In this paper, a short summary will first be given 
of the theory, which is then applied to the calculation 
of a practical distillation. 

It is impossible to present even an abridged version 
of the theory in the space available here. Therefore 
only the course of the reasoning is indicated and all 
equations that are of fundamental importance, or are 
used in the practical application, are given without 
proof. The reader, who is interested in the derivations, 
is referred to the original publication. 
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GENERAL EQuatTions 
Let us consider the stripping section of a distillation 
column. The trays are numbered upwards, so tray 
number m lies below the (m + 1)th tray. The reboiler 
is indicated by m = 0. 
The liquid and vapour rates of any arbitrary 


component are related by the basic equations 


, 9 
i,m iw ae x, at (<) 


From these equations one can calculate L, ,, , ; and 
doing this one obtains the well known expression: 


kom 
\ 1 I I | 
mil B.S 1 (A,,%, lien 1% m il , a me x) 3 
k-0 Pe 
4 1 1 
+ U, MA. a +» 4A, @; »)- 


The product of the absorption factors of the reference 


component on the successive trays 1,2... & will be 
written as Ap(k). So, 
Ap(k) = A, . A, ... A, 
and by definition 
Ap() =1 (4) 


Asimilar notation is adopted for the product of the 
relative volatilities: 
ap(r.O)—1. (5) 


apli, k)=a;,.a%9..-%, and 


Introducing the newly defined functions Ap(h) 


kom 
and ap(i, &) in eq. (3) and abbreviating Y’ to 2, leads 
k-0 
to: 
; l | yy aplt, m) ; , | : 
‘imti Ap(m) | Bee api, k) Ap(k) Vio %p(t,m) - (6) 


For a constant relative volatility one has «p(i, &) 
xt and then eq. (3) becomes: 


I { m | Loyd 


L ; B. p ae * Ap(k) I 0 x; j- (4) 


im 


Ap(m) 


A constant relative volatility will be assumed in the 
following. 

This equation expresses L,,,,, in terms of the 
absorption factor product for all components for which 
a residue rate B; and a vapour reflux V; » is specified. 
It is, however, not possible to specify beforehand 
both the distillate rate and the residue rate for all 
components in a multicomponent mixture. [It can be 
proven that this can be done for 4 components at 
most, so that there are a number of components for 
which either the distillate rate or the residue rate re- 
mains unspecified. Components for which both D,; as 


+t A list of notations is given at the end of the paper. 
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well as B; are specified are called key components. 
If the distillate rate is unspecified but diminishing as 
compared with the total distillate rate, the component 
is called unspecified in the top, and, if the residue rate 
is negligible, the component is unspecified in the 
bottom section. Components which have a residue- 
and distillate rate that cannot be specified @ priori 
but which are both non-negligible are called inter- 
mediate components. 

In the following, only two key components are 
assumed to be present and no intermediate components. 
There are, however, unspecified components. For 
those which have a negligible residue rate (they are 
volatile components) it is not practical to use eq. (6), 
especially since their residue rate varies with number 
of plates and reflux, though it remains extremely 
small. These components are present in finite quant- 
ities only in the top section and on a few trays below 
the feed tray. We shall, therefore, express their liquid 
rate in the absorption factor product and in their 
vapour rates at the feed tray. One obtains from 
eq. (1) and (2) if B; becomes equal to zero: 


o™ (M1) 


l . = 


hie ‘tym APM +) leu: @ 
The index M + 1 is used to denote the feed plate when 
it is considered as a tray in the bottom section and 
N +1 is used when it is considered as belonging to 
1) is the 
product of absorption factors in the stripping section 
on tray M+ 1, 

It has been shown that a distillation in which the 


the top section. The magnitude Ap(M 


which is the feed tray. 


reflux is variable can be calculated in exactly the same 
way as a distillation at constant reflux, if instead of 
the liquid and vapour rates other quantities are intro- 
duced. They are called generalized liquid and vapour 
rates and have been defined in the previous papers 
[5], (6). 


discussed here. 


For simplicity, constant reflux only is 


Constant Motar Reriux 

For a distillation in which the molar liquid rate may 
be considered as constant in the stripping section, 
addition of the liquid rates for all components leads to 
a constant value L,. The resulting equation is used 
for the calculation of the absorption factor product. 
One obtains on addition of eq.7 and 8 for all compo- 
nents : 

2; (2, Boat* Ap(k) 


b APM 1) pp Be a OOO, wa, 


Ap(m) L, Viom'} 


(9) 
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The notation Y, indicates that only the components 
with a specified residue rate are included in the sum- 
mation and 2” means a summation over the un- 
specified components. It may be recalled that besides 
constant reflux, constant ratio of volatilities is assumed 
and two keys and no intermediate components pre- 
One obtains for the solution of eq.(9) the 


following expression : 


Ap(M 4 


sent. 


Ap(m) = Ap,(m) 1) 2 Cat (M +1) (10) 


In this equation Ap,(m) has been written for the 
solution of eq.(9) if all unspecified components are 
omitted and (,, is: 

Vy VuNoa +. 


k—m 


; (11) 
ua Ly — = %y 


u(m—k) 
k) is written as an abbreviation of 


2 Bat, 


The function u(m 


At the feed tray the same temperature and, if the 
variation of the pressure along the column must be 
taken into account, the same pressure must be found, 
irrespective of the fact that its absorption factor is 
obtained by a calculation from the absorption factor 
product in the stripping section or from the top 
section. Therefore, the same equilibrium constant 
of the reference component must be obtained in both 
cases. This condition provides a relation between the 
numbers of ideal trays in both sections and the reflux 
ratio. The absorption factor in the stripping section 
at the mth tray is equal to A p(m) Ap(m_—1) and 
therefore, the reciprocal value of the equilibrium con- 
stant at the mth tray is equal to 
Ap(m) 


Pan Vy 
nm Ve Ap(m—1)~ 


(12) 
In the top section similar equations are obtained 
as for the stripping section if V;,,,,, L;,, Dj and &,, 
are substituted for Ly»... Vim, By; and A; respec- 
tively. The plates are numbered downwards and n = 0 
for the condensor. Instead of the absorption factor 
product, here the stripping factor product Sp(n) can 
be calculated from a linear summation equation, which 
is the analogon of eq.(9). One finds for a component 
which is specified in the top section: 
l 


Spin) { D, Z,az"** Sp(k) + Ly ga; "} (13) 


inl 

t In the first publication c, was used instead of C,. This 
might, however, cause confusion with constants c;, defined 
in that article, which have a slightly different meaning. The 
approximative value D, is used in that paper in the final 
equations, instead of V, y ,, which was used in the derivation 
itself [5}. 
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and for an unspecified component: 


LUN 
Vy _ Sp(N +1 Ly 7 
, . pt: rt) Lt “u,M+1° 


until 


(14) 


u 
Sp(n) 


The equation for the calculation of Sp(n) becomes: 
Sp(n) V, = 2; { y, Da ” Sp(k) 4 Li; 9 %& | 


. (V44) 7 ‘ re 
“u,! B* 


Sp(N +1) -s s Ce 

Of course the components which are unspecified 
in the top section are different from those in the 
bottom section and, therefore, a component with 
index wu in eq.(15) is different from the components 
denoted by uw in eq. (9). If it is necessary to point to 
this distinction, a second subscript 6 or t respectively, 
is added to indicate that the component is unspecified 
in the bottom or top section, respectively. 

The solution of eq. (15) is: 


Sp(n) = Sp,(n) + Sp(N +1) 2, 


for moderate and large values of N and 


lt Lu, M +1 
Ly V,— 3, a%* u(n—k) 


’ 
“ 


(17) 


k) is written for p > Da; *) Since 


1) is the stripping factor in the top 


Here u(n 
Sp(n) Sp(n 
section, the equilibrium constant of the reference 


component is obtained as: 


» Lt NSp(n) ' 
4 
Kron Vi Sp(n—1)° (18) 


At the feed tray we must have: 


Vp Ap(M +1) Ve Sp(N) 


= ( 
ly Ap(M) Lt Sp(N + 1° (19) 


It is clear from eq. (10) and (16), that first the ab- 
sorption factor product and the stripping factor pro- 
duct must be calculated, neglecting the unspecified 
components. This can be done either step by step, so 
that first Ap(1) is caleulated for m = 1, then this value 
is inserted in eq. (9) for m = 2 and Ap(2) is calculated 
and so on, or by solving eq.(10) analytically. The 
analytical solutions are: 


Ap, (m) = 2; Co Pir 20) 


Sp,(n) = Le; Pi e- (21) 


The second subscripts are added here in order to make 
it clear that different constants are obtained in both 
formulae, but they will not be retained in cases where 
confusion is out of question. The symbols have the 
following meaning: 
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The quantity g, } is a root of the algebraic equation 
containing the purely mathematical magnitude ¢ as 


unknown variable: 


B 
‘ ’ »» 
L=2 a (=2) 
lim Vio (lah) 
Chb (l— qo) ; 23) 
. t Vb , L, ~ B, (I a t) 
and g;; is a root of the algebraic equation 
, D 
‘ . » 
= 2, = 24) 
, 1 
| 2; L, (la) 7 
C: im , (I 4 ,t) 1,0 (25) 
(=9,, V¥,—-2: Dl — a; "t) 
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and 
Sp, (N) 
Sp, (N +1) 


%  Sp(N) 
la Sp(N , 1) 


iv | 


I—Z ¢ 
Lt ' 


=a ~S 
Discussion OF EXAMPLE 

The theory will now be applied to the actual calcula- 
tion of a distillation. A separation will be considered 
that has been treated by another author (JENNY) who 
used a plate-to-plate calculation method of the liquid 
and vapour rates [2]. The separation scheme is as 
follows: 

795°C. The 


relative volatilities %; » correspond to that temperature. 


The temperature of the feed is 


There is some uncertainty about the best choice of 
the average relative volatili- 








Table 1. Ne paration acheme of 6 comp. syste sadomt ssure 20 atm. The product rates are in ties in both sections owing to 
moles per unit of time — . 
the variation of «; with tem- 
Vapour Ti re ari tic “ 
Comp. Diet. Res. “a Sie. fond tas _ ma perature. The arithmetic mean 
' value of the top- and feed 
temperature has been taken 
l 250 O20 OOoLo 6 34- 
‘ pee as the relevant temperature 
» o-0o0 (O78 ool 5-07 7-2 teh é 
rf A for «; , and similarly the arith- 
3 247 0-004 180 O70 24K) 2.63 1-96 a ‘ 
- " ae etic 7 value twee 
4 0-003 0-167 0-095 O-OTS 100 1-00 1-00 metic mean 1 ret n 
5 O10 0-038 0-072 0-429 0-497 bottom- and feed temperature 
6 0-120 O-O25 OAS 0-206 0-272 has been used for Ai be This 
applies in both cases to the 
0-599 401 0-666 0-334 


specified components only. The 





The eq.(22) and (24) are, in principle, the same 
equations which have also been obtained by the authors 
referred to in the beginning of this paper. However, it 
follows from the structure of the eq. (10) and (16) that 
not the individual components but only the functions 
k) and u(n—k) the 


summation. This means that all distillations in which 


u(m appear under sign of 
these functions are the same, lead to the same eq. (22) 
and (24). In the authors’ opinion, this is a point of 
fundamental importance for the theory of distillation. 
It is also very important for practical application, 
since it opens the possibility of a well founded method 
for approximation of a multi-component mixture by 
& mixture containing less components. The degree 
of the eq.(22) and (24) may thereby appreciably be 
reduced [7]. 

After the functions Ap,(m) and Sp,(n) have been 
obtained a simple calculation shows that: 


Ve, Ap(M+1) Ve Ap,(M +1) 
ly <Ap(M) ly Ap iM) | 
" (26) 
1— 200, + Ap, (At +1) (LEC, a") Ap,(M) | 


relative volatility of an un- 
specified component will be taken equal to % » since 
these components occur only on a few trays near the 
feed tray. There exist theoretical reasons for a choice 
of the average « of the specified components which 
is closer to the relative volatility at the feed tray, 
if large numbers of trays are involved. 

In the following calculations the values from Table 1 
have, however, been used in all cases since it is thought 
that before refinements are introduced, the straight 
forward calculation methods should be discussed first. 
The functions u(m——k) for the bottom section and 
u(n—k) for the top section are represented in Fig. la 
and b. These functions have been obtained as the sum 
of the 4 straight lines which are the graphs of B,«" * 
or D,a; “"~” when plotted on semi-logarithmic paper. 
The dotted lines and the circles in the figures will be 


discussed later. 


STEPWISE CALCULATION OF Ap(m) AND Sp(n) 
The absorption factor product for any assumed reflux 
ratio can now be calculated from eq. (9) in which the 


unspecified components are neglected. For a reflux 
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ratio Rp = 1-50 in the top section one has L, = 0-899, 
V, = 1-498, L, = 1-233, V, = 0-832 and thus one 
has if V;9 is set equal to B,V,/B ic. if the vapour 
reflux has the same composition as the residue 
Ap, (0) = 1, Ap, (1) x 0-832 = 0-264 + 2-06 > 
0-264, Ap, (1) = 0-972. 

Ap, (2) x 0-832 = 0-972 > 
2-06 « 0-218, 

A p,(2) = L115. 

Ap, (3) x 0-832 = 1-115 x 0-264 +. 0-972 > 
0-218 + 0-213 + 2-06 x 0-213, 

Ap, (3) = 1-395 


u(m-k) 


0-264 + 0-218 4 to! 


and so on. For the calculation of Sp,(n) in 
the top section eq.(15) is used and we have 
Ly, 9 = D,;.1,/D im this case. The values of 
Ap,(m), Sp,(n) and V, Ap, (M+ 1) L, Ap,(M) 
and V,Sp,(N)/L,Sp,(N +1) are given in 
Table 2. The latter two variables have been 
denoted by §, and », respectively. They are 
equal to the reciprocal value of the equi- 
librium constant of the reference component. 
This is, however, no longer so with variable 
reflux. 

This stepwise calculation is different from 
conventional plate-to-plate calculations. In 
the latter material balance equations and 
equilibrium conditions are used alternatively 
for each component to calculate its flow rate or 
concentration. Here the single eq.(9) or (15) is used 
to calculate the absorption or stripping factor pro- 
duct. The equilibrium conditions are automatically 
satisfied. 

The unspecified components have been neglected 
in these calculations as is indicated by the subscript | 
added to Ap(m), Sp(n), § and y. 


ig. la. Characteristic function 
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CORRECTION FOR UNSPECIFIED COMPONENTS 
We shall now correct the above values for the presence 
of the unspecified components, assuming that the tray 
under consideration has been selected as the feed tray. 


10 


m-k n-k 


Fig. 1b. Characteristic function 
for top section. 


for stripping section. 
In that case § and » are given by the expressions 26 
and 27 respectively. The constants (', are calculated 
by eq. (11) for the stripping section and by eq. (17) for 
the top section. Owing to the fact that the relative 
volatility of an unspecified component is high in the 
bottom section and that it is small in the top section, 
important simplifications can be introduced in these 
formulae. 





Table 2. Absorption and stripping factors and §,n values for 6 component system at R,, = 1-50 


morn Ap, (m) Sp, (n) &,(m) H, (nm) 
0-656 
0-799 
0-844 
0-890 
0-939 
0-994 
1-047 
1-110 
1-179 
1-232 


0-328 

0-156 

0-0822 
0-0465 
0-0293 
0-0211 
0-O0171 
0-0151 
0-0142 
0-0137 


0-972 
1-115 
1-395 
1-840 
2-560 
3-770 
5-849 
9-687 
16-920 
30-880 


&,;(M +1) ny (N +1) &,,;(M +1) nN +1) 


5-080 
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If all terms containing «,' are neglected in eq. (26), 
one finds in the stripping section for the function & 
corrected to the first order: 

2” Dy 


E,(M +1) = §,(M+1)/(1 ; (28) 
t 


ve 


Here 2” D, 1-498 ; 


therefore £,(M 


dD, 
1) 


D, = 0-349 and |, 
1-30 &,(M + 1). 


fOR1 


5 


8 10 ° 
M+1 OR Nel 


Fig. 2. Values of = and » for six component system as 
a function of M+ 1 and N+ 1. 
lation for Rp = 1-50. 
II. corrected to second order. 


Stepwise calcu 
1. not corrected for unspec. comp. 
Analytical calcula 


tion, not corrected for unspec. comp. 


| 
| 


Analytical calculation, corrected to second order. 


1-04 


78 
12 


765 


e Rp 1: 
tla Rpt 
v Rp 


Ry 
Rp -~ 
Rp 0-764 


Ry 0 51 


In the top section, the first order correction of 4 is 
obtained by neglecting in eq. (27) all terms that contain 
a factor «,. This gives: 

2” B 

ly 

B; + Be — 0230 and L, 
1) = 0-744, (N +1). 


The values of & and » corrected to the second order 
1 


ny (N+ 1) "\ny(N+1). (20) 


vv 


Since 2” B, 1-899, 


one has 4,(N 4 
are obtained when terms containing a factor a * in 
eq. (11) and a factor «, in eq. (16), respectively, are also 
retained but higher powers are still neglected. One 
has: 


E,,(M +1) 


1) 

(30) 

, 2 Bain 
Vy 


Ap,(M +1) 
Ap,(M) 
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nil ; 1) my (X t 1) » 


2” B, uF 
Le 


"Dip \) 
le II 
Vt 
“In j 
xy" —> 0-030 , 
0-0726, 2” D, asp 


(31 
=” By Gu, F 
Lo 
In the 
=” B; a; , 
0-153. 
and (31) one finds 


5), (M +1) =&(M 


present example 
ont. _ 
0-247 and 2 By a, » 


If these values are inserted in eq. (30) 


1) 0-020 0-600 & 
Ap,(M +1) | 


ApiM) |} 


767 

| 
0-206 

and 

1) 


1) (0-744 O-050 B75 . WmI7TL) & 0-099) ‘ 
J 


In these corrections it has been assumed that the 


Hy N 
Wy (XN ; 


product rates of the unspecified components are 
negligible compared with those of the specified compo- 
nents. This condition can only be fulfilled by a suffi- 
cient number of trays in each section. It is for that 
reason that only corrected values of & and 4 from m 
or # <= 4 upwards have been included in Table 2. The 
amount of unspecified components in the terminal 
products can be estimated from eq.(8) or (14). 
The values § and » have been plotted as a function 
of M+1 or N 
(Fig.2). The feed may be introduced at any tray 
for 


| respectively in the same graph 
which § =. This graph gives therefore all 
columns with which the separation can be performed 
for the selected value of L, if the calculation of Ap, (m) 
and Sp,(n) is continued sufficiently far. Possible 
eg. M+1=10 N+1=7 
or 9 trays below the feed tray and 6 trays above it, 
or M+1=N+4+1= 85. 


combinations are and 


Analytical solution 
The functions Ap, (m) and Sp,(n) may also be caleul- 
ated from the analytical solutions eq.(20) and (21) 
respectively. The eq.(22) and (24) are both of the 
fourth degree in our case and solving these would be 
a laborious procedure. We shall therefore approxim- 
k) by those for 


m 


6 
B, «3 + B, is approximated by a single exponential 


ate the functions u(m-——k) and u(n 


a binary mixture. For this purpose the sum B, « 


function B, a" in the stripping section from m = | 
on upwards [3]. The latter function is represented 
by the dotted straight line in Fig. 1a; one has «, = 0-88 
and B, = 0-253. The function u*(m) = B,a}" + Byaf' 
is now an approximation for u(m) from m = 1 up- 


wards. Introducing L} = L, + u*(0)—u(0) and 
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Table 3. Data for analytical solution of eq.(9) 





» , 
L} Ly V;, le 


1078 
0-666 
0-456 


1-402 
1-000 
0-790 


1-001 
0-509 


1-258 
O-S56 


0-646 O-389 


1-78 
1-12 
0-765 


Rp P3,b Pb C3.b Ch 


00-0219 
0-0375 


O-O775 


1-97 0-69 
1-98 


2-00 


1-10 
1-25 


1-43 


0-71 





inserting u*(m-—k) instead of u(m—k) for m —k 


1, 2, 3 ete. in eq. (9) leads to the equation: 


Ap,(m) LF = XS, u* (m R,,u* (m). (32) 


k) Ap, (k) 
The unspecified components are again neglected 


m 
Cy Pr 


m==1 on upwards, whereas Ap,(0) = 1 


here. Its solution is Ap,(m) + €3q 3 from 
and q,' 


and q,' are the roots of the quadratic equation 


* I 3h 
*b l1— apt 


The constants ¢ are according to eq.(23) equal to 


Ly 39>) 


u(O)~ l 


t, Py y( 
1 
3th 


| Ag (1 


u(O) l 


a, ps *) (I 3% *) 


Vn%s3 . 

Instead of assuming a value for L* and solving 
the quadratic equation it is more convenient to assume 
a value for one of the roots say g,* and calculate LF 
and g,* which can be done by solving a linear equation 
in either case. One has gy, ' gy,’ = (L — B, — B,)/agx%. 
The data for these calculations for some L* values are 
given in Table 3. 

The following values of Ap? (m) have been cal- 
culated from this Table. 


Table 4. Values of Ap (m) 





Lt Apt(4)  Ap*(5)  Ap*(8)  Aps(9) — Ap#(10) 


21-08 
41-30 
100-00 


11-23 
22-80 


57-40 


6-47 
13-03 
32-20 


1-343 
2-331 
4-200 


1-76 
3°34 


6-72 


1-258 
0-856 
0-646 





The corresponding values of £, are plotted in Fig. 2 
and also the values of £,, corrected to the second order. 


~~. 
Ly 73, b- 


It follows from eq.(20) that §&,(m = oo) 
In the top section u(n) is approximated by 

u* (n) Da," 

a, = 2-76 and D, 

u(V) 

Sp, (n) 


Sp, (n) V,* =2,u* (n—k) Sp, (k) +Rpu*(n) n=1,2,3... 


D, from n =1 on upwards. Here 
0-287. Introducing y* V, 


u*(0) one obtains the following equation for 


C, ve Cy Ys and Sp, (O)—1. 


Its solution is Sp, (n) 
The equation 
y* ay) dD, 
é =A, " 1-8 
l—a,"t 
has the roots g,' and g,' and the constants ¢, and ¢, 


are equal to 


and 


-. 
1 % CI 


u(0) l %eVa I 

The roots q, and g are related by the expression 
G,'% | ~ Vi —D,— Dy ay. The calculated values 
of g, ete. are collected in Table 5 and the 4's are shown 
in Fig. 2. 

In graph 3 the top reflux Rp) has been chosen as 
the independent variable. The curves are — and 
values for a constant number of trays in each section. 
The corrected curves for an infinite number of trays 
intersect at Rp = 1-0. This is the minimum reflux 
for which Jenny has found 0-96 by plate to plate 
calculations. Differences of a few percent may be 
expected owing to the special choice of the relative 


Table 5. Analytical solution for Sp¥ (n) 





Ly 


1-160 
0-760 
0-457 
0-305 


Fat Ft Cat Ct 


0-575 
0-615 
0-660 
0-670 


0-01023 
0-0144 
0-027 
0-060 


0-450 
0 497 
0-585 
0-690 


1-003 
1-005 
1-010 
1-020 
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M+1=9 


1-5 so 


volatilities. The curve intersects the 
curve N+1=9 at Rp that 8+ 148 
trays are necessary to perform the separation at this 
1475 


reflux ratio, whereas we have found 7-5-4 


by the stepwise solution method. Graph 3 gives a 





OR" 
a 
i 
+ =~ 
a es 
: 1 
: 
o 
—_ 
2 
—— 
=8 
. 
1 
0S 2 
Rr 
— 
Fig. 3. Values of £ and » for six component as a function 
of Rp. — and » not corrected for unspec. comp. 


from analytical calculation. 


4M+1 5 oO N l 5 
V7M+1=9 oN l 9 
4M+l=10 @N-+1!1 10 
V7M+1 x ON l x 
wees idem, but corrected for unspec. comp. to second 
order. . € and » not corrected for unspec. comp. 
from McCape-Turece diagrams. 
aM+1=5 oN+1 5 
17 M+1=+9 ON ] a 
“leo M+1=00 NeNn+l=10 
eW-1 x @N+1 x 


and =» = for M+ 1 7; 
N+ 1 7. 


+» for N+ 1 8; 


complete survey of all possible columns and all the 


reflux ratios with which the desired distillate and 
residue can be obtained. 

Another useful representation is given in Fig.4 in 
which the curves represent € and 7 values for a constant 
value of J, and J, respectively, and the number of 
trays in each section is plotted as absissa. This re- 
presentation is still more general than the former one 
since it can also be used if the degree of vaporization 
of the feed is varied. This simply means that another 
I, corresponds to a given V, and so one has to 


combine different pairs of curves for different degrees 
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of vaporization. So for example a minimum reflux 
of Rp = 0-67 would have been found in case of no 
vaporization of the feed Vj, = V,. The representation 
in Fig.3 could also have been given a more general 
character if instead of Rp the corresponding J, and V, 


had been plotted. 


EOR* 
5 
4 P<. 
» : i. 
+ ae es aa | 
es ee 
z- a = a ~~ 
| a = ~~ . 
1 . ——— ve 3 , 
% 2 4 6 & io~ 
M+1 OR Ne! 


Fig. 4. Values of £ and » for six component system as 
a function of M+ 1 and N + 1. Not corrected 
for unspec. comp. from analytical calculation. 


@ L, = 1-402 o V, = 1-759 
E40 L, = 1-000 OV, = 1-359 
v L, = 0-790 16 V, = 1-056 
4 V, = 0-904 


eee &and 7» corrected to second order for unspec. comp. 
- € and 7» not corrected from McCanr-THrece dia- 
grams. 


GRAPHICAL SOLUTION 

The simplified eq. (32) is very similar to the equa- 
tion for the binary mixture of the components 3 and 5. 
The only difference is that Ry appears in eq. (32) in- 
stead of 

Ry, = {LF —u* (0)} /u* (0). 

It can be shown that the absorption factor products 
calculated from eq.(32) are (R, +- 1)(RE +-1) times 
the absorption factor products of the binary mixt- 
ure [8]. Therefore, the absorption factors of the latter 
are equal to Apt (m+ 1)/Apf (m) if m exceeds zero 
and the absorption factor at the first tray in the 
binary mixture is equal to 


{(RB +1)/(Rg +1)} Apt (1) 


since Ap(0)=1 in each case. The absorption factors 
of the binary mixture can be found easily by one of 
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the current graphical methods, from McCaBe and 
THIELEs diagram, for instance. This diagram is given 
in Fig. 5a for the bottom section and in Fig.5b for the 
top section. The equilibrium line in Fig.5a has been 
drawn for the components 3 and 6 which have a rela- 
tive volatility of 1-96/0-88 = 2-23. The concentration 
of the volatile component in the residue is 2, 

0-004 /{0-25% 0-004} = 0-0156. So the operating 


wo 


: eee we a | 
40 60 60 100 
Jo LIQUID OF COMP 3 
Fig. 5a. Separation of components > and 3 in stripping 
section, relative volatility = 2-23. 
for Rp 1: 
Rp L-1: 
Rp 7 
Indication of the trays 
® belonging to Rp 
Vv belonging to Rp 
® belonging to Rp 


0-0156 and x 


The ordinary 


lines are drawn across the point y 
0-0156 and their slopes are L*/V;* 
construction is then carried out to arrive at the con- 
centrations at any arbitrary tray. If x,, and y,, are 
the concentrations of component 3 in liquid and vapour 
respectively leaving the mth tray, one has Ky}, 

Tm/Ym at that tray and the absorption factor for 
component 3 is L*/V,* K,,,. This absorption factor 
is equal to that of the same component in the simplified 
multicomponent system. Therefore, one obtains &,(m) 
for this latter system, by multiplication of this ab- 


sorption factor by V,/L,a;'. So we have: 


(m) = (V, LE/L, Vi") +g <(a/Y%_) m= 2,3,... (33) 


= 
- 
> 


The value of K,}, is read at the upper horizontal 
line in the figure as the intersection with a straight 
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line through the origin and the point representing 
the concentrations at the mth tray. So, for instance 
one reads for L* = 1-258 K,} = 0-704 and therefore 
Apt (9)/Apf (8) = 1-73 and &,(9) = 1-235. This tallies 
with the value 1-74 calculated from Table 5. As can 
be seen from Figures 3 and 4, all € and 7 values caleu- 
lated from the McCape-TuHreLe diagrams agree ex- 
cellently with the other values. 
100 ; - 


, 


f 


MOL io VAPOR OF COMP t 


MOL %/o LIQUID OF COMP L 


Fig. 5b. Separation of component ¢ and 4 in top section, 


relative volatility = 2-76. 
for Rp = 1-94 
for Rp 1-27 
for Rp = 0-764 


for Rp 


Indication of the trays 


@ belonging to Rp 


@ belonging 


VY belonging 


to Rp 
to Rp 


1-04 
1-27 
0-764 


Y belonging to Rp 


{pr* 
: V 


u* (O)} /u* (0) 
and the stripping factor product of the binary mixture 
. * 
is (Rp 
the multicomponent mixture. The stripping factors 


In the top section one has RF, 
1)(Rp, +1) times the stripping factor of 
Sp(n-+-1)/Sp(n) are therefore equal in both cases 
except for n = 0. 

A McCase and Tarece diagram has been con- 
structed for the components ¢ and 4 (Fig.5b). Their 
relative volatility is o/x,— 2-76 and y = 0-287 
(0-287 +. 0-003) = 0-990. Operating lines have been 
drawn for the vapor rates in Table 4. Their slopes 
are L*/V,*. The absorption factor of component ¢ at 

al *)77* 7 
the nth tray is L//V,"K, ,, 
absorption factor of the reference component 4 is, 


and K,;) = z,/y,. The 


therefore, ia ' 
(LF /V,") (Xn! Yn) % 
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This gives for », (»): 


(34) 


(V, L/L, Ve") ay (ey/y,) 0 


1, (m) 


From Fig.5b one reads e.g. for V,* = 0-747 at the 
Mh tray 2,/y, ~ 0647. This gives »,(9) — 2-51, 
which again tallies excellently with the analytical 
calculation. 

A list of and » values calculated by this graphical 
method is given in Table 6. 


Table 6. §&, and », values from McCane-Tutere diagrams 





Rp=112  Rp= 0-765 
951 
0-898 


872 


1-15 
1-06 
1-03 
0-904 
0-951 
0-900 
0-855 
0-814 


O-S56 
O-S834 
0-813 
0-787 
0-761 








127 Rp =0-764 Rp=0 


2-84 
315 
3-26 
3-39 
3-4 
3-68 
3-82 
3-04 


2-24 
2-40 
2-52 
2-68 
2-90 
3-14 
3-37 
3-56 
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NOTATION 


The standard system of nomenclature for chemical engineering 
operations (Trans. Amer. Inst. Chem. Engrs 1944 40 251) has 
been followed as far as possible. A* indicates that the quantity 
refers to the binary mixture approximating a multicomponent 
mixture. 

x, (alpha) — Relative volatility of component i towards a 


reference component. 
k= 

A= ba Ap(k) Generating function of Ap(m). 
k=0 
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» Absorption factor for component i on tray m. 
Absorption factor for reference component 
on tray m. 

A, A, ...A,,- Absorption factor product of 
reference component on tray m. 

Residue rate, moles/sec. 

Residue rate of component i, moles sec. 
Subscript referring to bottom section. 
Constants in analytical expression of A. 
Constant used in correction for unspecified 
components, 

Distillate rate, moles/sec. 

Distillate rate of component i, moles/sec. 
Root of equation to be solved to obtain poles 
of A or 8. 

Index denoting a certain component. 
Equilibrium constant for component i at 
tray m. 

Liquid rate below feed, moles/sec. 

Liquid rate above feed, moles/sec. 

Liquid rate of component i from tray m, 
moles sec. 

Index for tray below feed, m = 0 means 
reboiler. 

Index for feed tray. 

Index for tray above feed, n = 0 means 

condensor. If in the symbols m is replaced 

by n they refer to the top section. 

Number of components. 

Index for feed tray. 

Reflux ratio in top and bottom resp. 
ko 
> Sp(k) ® Generating function of Sp(n). 

k=O 

Stripping factor of comp. i at tray n. 
Stripping factor for reference component on 

tray n. 

S, S,... 8, Stripping factor product of re- 

ference component on tray n. 

Variable in generating functions or subscript 


referring to top section. 


(m hy) 
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k=co 

S u(k)t*. Generating function of u(m—k) 

k-0 

or u(n—k) respectively. 
Vapor rate below feed, moles/sec. 
Vapor rate above feed, moles/sec. 
Vapor rate of component i from tray m, 
cal /sec. 

ims Xp = Molal liquid concentration of component i at 

tray m or in residue, respectively. 

E (xi) = Ap(M +1) Vy)/Ap(M) Ly), parameter for 

feed tray. 

Yim: Yip = Molal vapor concentration of component i 

at tray m or in distillate, respectively. 

n (eta) — Sp(N) V,/Sp(N+1) Ly), parameter for 


feed tray. 


p= % 


m .« . 
2" in bottom section. 


n: : 
; In top section. 
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Summary 


The free energies of formation (.1@) of the fluorides, chlorides and oxides of the alkali, alkaline 


earth and rare earth metals are inferred tentatively from available data, and the results are given in the 
form of “‘Ellingham”’ diagrams ** in which 4G is plotted against 7’. 


From the AG values the relative tendencies of the various chlorides to convert a fluoride or oxide to 
the corresponding chloride are deduced and the results found to accord with experiment in cases for which 


information is available. 


Indications have been given elsewhere [2], [3] of the 
interest to the magnesium industry of a knowledge of 
the positions of magnesium in the chloride and fluoride 
electropotential series, since these positions will largely 
govern the choice of possible constituents of salt 
mixtures for alloying and purifying the metal. The 
order of metals in the two series is different, the dif- 
ferences arising from variations in the magnitude and 
sign of the heat of solution of the salts concerned. In 
the present paper an attempt is made to infer the free 
energy - temperature relationships for the chlorides, 
fluorides and oxides of the alkali, alkaline earth, rare 
earth and certain other metals from published data 
on decomposition potentials and exchange reactions, 
supplemented by approximate thermochemical calcu- 


* This paper is based on work forming an appendix to 
a thesis approved by the University of London for the degree 
of Ph. D. (Thesis submitted Ist May 1950.) 


** This convenient term has been used by Parsons [1 }. 


lation. Consideration is also given to the occurrence 
of double decomposition reactions involving these 
substances. 

Since drafting this paper, Ellingham diagrams 
based on calculations of the extended type have been 
published for chlorides by OsBorn [43], KELLOGG [44] 
and Vina [54], and for fluorides by KeLLtoee [55}. 
Data for chlorides 
published by Brewer, Bromiey, Gittes and Lor. 
GREN [45]. 


and fluorides have also been 


OsBorN shows curves for four of the chlorides of 
present interest, namely, NaCl, MgCl,, CaCl,, and 
BaCl,. The relative positions given for the first three 
do not accord in any way with the conclusions of the 
present papert. On the other hand the results of 


+ Note added in proof: It has since been explained by 
OsporN [78] that 50 keal should be added to the curves for 
NaCl and CaCl,, and this evidently applies also to that for 
BaCl,. Osporns results have been added to Table 10 
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KELLOGG and Vita for these chlorides, and for the 
other halides of present interest with which they deal, 
are in reasonably good agreement with the results 
here deduced (Table 10). 


Brewer, Bromiey, Gittes and LorGren alone 
give data for all the chlorides and fluorides of Tables 2 
and 4. Their figures have been selected critically, but 
apparently without consideration having been given 
to the available experimental data on decomposition 
potentials. These figures are not therefore thought 
likely to give a truer picture of the relative stabilities 
of the chlorides and fluorides involved than the data 
deduced in the present paper, particularly since in 
the case of the fluorides, where decomposition poten- 
tial measurements were not available for guidance, the 
more exact calculations of KeLLoGe are in rather 
better agreement with the rough calculations of Table 2 
than with Brewers figures. A number of entropy 


determinations havem oreover become available since 


Chemical 
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the latter were deduced. In one or two cases the 
BREWER entropy estimates appear to be somewhat 


in error’. 


Tue FLvorRiIpE ELECTROPOTENTIAL SERIES 
Published information 


No satisfactory measurements of the decomposition 
potentials of fluorides appear to be available. Campi 
and Devoro [4] were unable to achieve electrode 
reversibility; and the earlier measures of NEUMANN 
and Ricurer [5] were made by an unsatisfactory 
method [6]-[8]. The results obtained by Neumann 
and RicuTer are numerically below the corresponding 
figures found for the chlorides and evidently much too 
low, despite agreement between the calculated values 


-AH) and 


of the heat of formation per mol (Q, 


* Topps [62] recent experimental determination of 13-68 
for MgF, is considerably below Brewers estimate of 18, 
though it agrees with a published estimate [50] of 13-6 based 
on semi-empirical methods [56], [57]. 


Table 1, Calculated approximate free energy of formation of various fluorides 





AH™** 
Substance (— Q,) 
(keals) 


v208 
Sur, * 


(E.U.) 


8-50) 


12-5) 


145-57 
135-95 
134-51 
263-8 
290-2 
289-0 
287-9 

4 ygi@a) 
320-000 
399-6" 


SSL? 


15-8 

13-6") 

16-4 

20-6"? 13: 
23-03 15- 
10) 2.8 
125" 6-73 
35-4! b) 


13-8 
ma) 13-6" 


BaF, 
BeF, 
AIF, 
CeF, 
Thr, 


1 10298 (Free energy 
n per qm. equiv) 


; (AH** — TAS**) 
(keala) 


A Ns 


+208 
SMP, 


(E.U.) 


138-9 
128-8 
127-5 
125-5 
138-8 
138-4 
137-9 
1143 
103-0 
128-1 
133-6 





* Unless otherwise stated these figures have been taken from Brcnhowsk1 and Rosstn1 [17]. 

+t Unless otherwise stated these figures have been taken from WENNER [56]. 

‘*) Calculated from data of TreapweLt and Mauper i [57] and in close agreement with the values given by Drosspacu 
(op. cit. p. 24) and Giasstone [58]. (A subsequent determination of 15-91 E.U. has been made for KF by Westrum and Prr- 


ZER [59)). 
(b) 

WENNER [56) and subsequent literature. 
‘) Estimated by Gross [60]. 


(4) Given as —332-0 to —500-0 by Renavutt (61). 


Estimated from plot of S against log. mol. wt. for all metal halides of the same valency, based on data given by 


(©) The figure given is the sum of the heats of formation of the aqueous ions. 


() Estimated from Latrmers Formula [63}. 


‘®) Estimated by method () and also by application of the method of Treapwett and Mauperwi [57] using the known 


values for CaF, and BaF,. 
(1) Metals Reference Book [64]. 
® Crustus, GoLpMANN and Perwick [65]. 
Note: The sign 


means here “numerically less than.”’ 


(A subsequent determination of 13-68 E.U. for MgF, has been made by Topp [62)). 
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calorimetric determinations. This agreement may be 
attributed to the very high numerical values found for 
the temperature coefficient* of the decomposition po- 
tential, dE/d7’, the effect of which has been to offset 
the low figures obtained for the decomposition po- 
tential** itself. AGar and Bowpen [9] showed that 
there are three main causes of irreversibility in elec- 
trode reactions: ‘‘concentration potential,’ “‘activa- 
tion overpotential’” and “resistance overpotential.” 
Since these quantities decrease with increasing tem- 
perature and dE dT is negative, all three sources of 
irreversibility would tend to increase the numerical 
value of dE/dT, and the high values found by Nev- 
MANN and RicuHTER are probably attributable to irre- 
versibility. 

Exchange experiments between magnesium and 
mixtures of alkali and alkaline earth metal fluorides 
have been made by Rurr and Buscu [11], and from 
these they infer that the electropotentials of the me- 
tals in their fluorides decrease in the following order: 
Ca—Li—Ba—Mg—Na—kK, thus representing a rever- 
sal of the usual order so far as Mg, Na and K are con- 
cerned. This inference, however, involves the assump- 
tion that the alkali metals are associated in molten 
magnesium, as was found to be the case in molten tin, 
antimony and bismuth [12]: otherwise the order 
Li—Ca—Ba—. Other difficulties in 
interpreting such experimental results are the lack of 


would begin 


information on activity coefficients or on the free 


energy of mixing or of compound formation in the 


metal phase. 

Although DrossBacu [13] considered the decompo- 
sition potential of NaF to be about 0-7 v above that 
of MgF,, there is no doubt that NaF and KF are read- 
ily reduced by molten magnesium; if these fluorides 
are also molten the reduction is violent with NaF and 
mildly explosive with KF. The necessity of removing 
NaF from all-fluoride welding fluxes for magnesium so 
as to avoid sodium introduction into the metal [14] 
and of the need for special steps to avoid explosive 
reaction between K,ZrF, and magnesium [15] are 
dE/dT = ASinF, where AS is the entropy 
increase accompanying formation of | mol of halide of valency 
n and F is the Faraday, and since the AS values for the forma- 
tion of the fluorides and chlorides are similar (see below), 
the values of dE/dT found for chlorides and fluorides should 
not differ widely; and yet NEUMANN: values for fluorides are 
many times those obtained for the chlorides, his figures for 
MgF, for example, being seven times that found by Cowen [10] 
for MgCl,. The latter’s figures are in good accord with calcula- 
tion. 

*Q,= E+ dE. 


* Since 
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practical reflections of the reducibility of these fluor- 
ides. The reduction of KF by magnesium has more- 
over been suggested as a method for the production of 
potassium f. 

The position of aluminium in the fluoride series 
appears to require clarification [2]. 


Thermochemical calculations 
Before attempting to infer the position of magnesium 
in the fluoride electropotential series, it may be of 
interest to consider the matter briefly from the 


Temperature ——= *K 


— 900 7100 1300 


“600 “B00 7000 
Temperature - a 





Fig. 1. Ellingham diagram for fluorides based on rough 
thermochemical calculations. 


Note. The Ellingham diagrams in Figs. 1-5 show the tem- 
perature variation of the free energy of formation of the 
halide or oxide from one gram equivalent of metal. 


thermodynamic point of view. Values of the free 
energy of formation per gram equivalent calculated 
from experimental or estimated data are given in 
Table 1. 

The AG T curves for reactions such as the forma- 
tion of the metal halides from their elements are ap- 
proximately linear above room temperature, and 


a linear equation will frequently represent the 


+ The low boiling point of potassium apparently enables 
the metal to be produced by reduction of KF even with 
aluminium [16)}. 
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observed data with the accuracy of the more comp- 
licated form of expression obtained by integrating 
the equation showing the temperature dependence 
of AC,. The justifiability of straight AG/T curves 
in the higher temperature ranges is illustrated by the 
linear dependence on temperature found by all workers 
for the decomposition potentials of various chlorides. 
If the necessary high temperature data were available 
in all cases, the AG/T relationship could be computed 
in the usual way and the results then represented by 
a linear equation for convenience. But such data are 
available only for some of the fluorides in the table 
and the values of AG,,, for even these are in some 
cases so uncertain* as to render attempts to apply the 


* The calculations have been revised to incorporate recent 
entropy determinations where available. 
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high temperature data of doubtful value. This is 
particularly true of the present investigation where 
the differences in the free energies of formation of the 
various fluorides are of more concern than the absolute 
values involved. A uniform procedure applicable 
to all cases seems therefore desirable, and for the 
present purpose it will suffice to indicate the tempera- 
ture dependence of AG by straight lines of slope A Sy, 
with inflections corresponding to the entropy changes 
accompanying changes of phase of the fluoride or its 
constituent elements. The results of such calculations 
are indicated in Table 2 and Fig. 1. 


Probable electropotential order of the fluorides 


It will be seen that, except for lithium, the calculated 
order of electropositiveness of the metals in their 


Table 2. Calculated approximate temperature variation of free energy of formation of various fluorides 





1H 1S Phase transformations AG/T equation . AG = A + B(T—C) 
Substance Table A Table A ‘ * 
“ga Temp. Nature of AStrars Range of T. I 
: c.U. Ms : " J } Ache 
(boat) (E ( RK) transform. (= L/T) ( K) A 1000 B ‘ n ie 
LiF 145-57 22-5 298— 459 145-57 22-5 0 
459 Li melts 1-66 459-1133 135-24 24-16 459 
1133 LiF melts 525 1133-1300 118-95 18-91 1133 115-8 
NaF 135-95 24-0 298— 371 135-95 24-0 0 
371 Na melts 1-70 371-1151 127-05 25-7 371 
1151 Na boils 20-2) 1151-1263 107-0 45-9 1151 
1263 NaF melts 6-18 1263-1300 101-8 39-72 1263 100-4 
KF 134-51 23-7 298— 337 134-5 23-7 0 
337 K melts 1-92 337-1031 126-53 25-62 337 
1031 K_ boils 29”) 1031-1119 108-75 45-62 1031 
1119 KF melts 5-56 1119-1300 104-73 40-06 1119 97-5 
MeF, 263-8 42-76 208- 924 263-8 42-8 0 
924 Mg melts 1-95 924-1300 224-25 «944-75 924 103-7 
CaF, 290-2 42-1 208-1124 290-2 42-1 0 
1124 Ca melts Psa — 1124-1300 242-88 43-04 1124 117-8 
SrF, 289-0) 41-3 298. 1044 289-0 41-3 0 
1044 Sr melts 1-37 1044-1300 245-9 42-6 1044 117-5 
BaF, 287-9 40-67 298 977 287-9 40-67 0 
977 Ba melts 1-51 977-1300 248-17 42-18 977 117-3 
AIF, 329-0 67-1 298— 933 329-0 67:1 0 
933 Al melts 2-47 933-1300 266-4 69-57 933 80-3 





* While the 4S transf. values have been calculated using 


heats of transformation and temperatures given by Bicnowsk1 


and Rosstn1 (17), the transformation temperatures given have mostly been taken from West [66] and Kaye and Lary [67 }. 


() 8. GLassTone, op. cit., p. 449. 
(>) The figure is assumed to be the same as for sodium. 


(©) Metals Reference Book (64). The figure of 8-95 E.U. corresponding to the latent heat of fusion given in references [68 )} 
and [17] is almost certainly incorrect, being not only widely different from the entropies of fusion of magnesium, strontium and 
barium, but lying outside the normal range of values found for elements. 
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fluorides agrees above 850° C with that inferred for 
800°C from the exchange experiments of Rurr and 
Buscu. Such a difference in the position of lithium 
in the two series is easely explained, since the curves 
for lithium, calcium, strontium and barium evidently 
lie ciose together, and lithium differs from the other 
three metals in forming solid solutions with mag- 
nesium and not compounds. When molten mag- 
nesium and NaF at temperatures above about 750°C 
are mixed in open crucibles, heat is evolved resulting 
in extensive local boiling of sodium. The percentage 
reduction reached under such conditions must there- 
fore considerably exceed the true equilibrium per- 
centage at 750° C, both because of the loss of sodium 
by volatilisation and also because the heat of reaction 
will have raised the temperature to a higher value at 
which a greater percentage reduction would be 
obtained at equilibrium. For practical purposes then, 
the result of the heat of reaction will be equivalent to 
lowering the temperature at which the kink in the 
free energy/temperature curve occurs by some 100°C 
or more, according to the scale of the experiment ; and 
such an effect may even be present to some extent in 
closed bomb experiments such as those of Rurr and 
Buscu. Similar remarks apply to the KF curve. The 
rough calculations indicate therefore that the fluoride 


electropotential series inferred from Rurr and Buscus 


results is in general accord with expectation, and it 
is not necessary to assume that it results primarily 
from complex formation as they suggest. 

Table | shows that the cerium rare earth metals 
evidently lie very close to magnesium in the fluoride 
electropotential series, and experiments have in fact 
confirmed that the rare earth metals are slightly more 
electropositive at normal alloying temperature*. The 
same applies to thorium. No information appears to be 
available concerning the position of uranium in the 
fluoride series**. Common membership of the actinide 
series might suggest a similar position to thorium, and 
heat data [17], [76] indicate that this is the case in 
the chloride electropotential series. 


ThCl, : }4H3,, = — 71-1 keal 
UCI, : 44H}, = — 71 keal. 


Insufficient data appear to be available to estimate 
the free energy of formation of ZrF,, but it is evident 


* This conclusion is based on exchange experiments 
assuming ideal behaviour. 

** Recent work by SAVERWALD [18] indicates that uranium 
resembles chromium, molybdenum, and tungsten in being 
insoluble in magnesium. 
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that the curve for ZrF, would lie well above that 
for MgF, on Fig. 1 since ZrF, can be reduced by 
sodium [19]. Moreover no unreacted ZrF, has been 
detected in magnesium-zirconium alloys made by the 
“inert fluoride’ process [20]. 


II. THe Caioripe ELecrroporentiaL SERIES 
Published information 

For the present purposes, a knowledge of the relative 
decomposition potentials of the various chlorides is of 
more importance than the establishment of the true 
reversible values for a few of them; consequently in 
considering the available information, attention will 
be directed mainly to the work of those who have 
determined decomposition potentials for a number of 
the chlorides which are of interest, even where it is 
probable that systematic errors are present and more 
accurate determinations on particular single salts 
have been made by others. 

The AGT curves for the various chlorides cal- 
culated from the decomposition potential data of Nrv- 
MANN and RicuTer [21], and Campi and Devoro [22] 
are shown in Figs. 2 and 3+. NEUMANN and RIcHTERs 
figures are markedly lower than those of Cams and 
Devoro; this is to be expected for the reasons already 
indicated in discussing the information on the fluoride 
electropotential series. On the other hand their 
values for dE/dT do not appear to be improbably 
high as was the case with the fluoride determinations, 
and it seems likely that the low decomposition po- 
tentials are here attributable mainly to the inherent 
defects of their method rather than to highly irre- 
versible electrode behaviour. In these circumstances 
it is probably fair to include NEUMANN and RICHTERs 
measurements when considering the relative electro- 
potentials of the chlorides but to exclude them when 
questions of absolute magnitude arise. Some equilibria 
in metal/metallic chloride exchange reactions involv- 
ing alkali and alkaline earth metals have been measured 
{23]-[27], but in the absence of information on ac- 
tivity coefficients and on compound formation, partic- 
ularly in the metal phase, there is no basis on which 
such results can be satisfactorily treated, and measure- 
ments of decomposition potential appear preferable 
for inferring the chloride electropotential series. 


Thermochemical calculations 
The results of thermochemical calculations of the free 
energy of formation of chlorides are given in Tables 3 


+ In converting the potential measurements to calories the 
Faraday has been taken as 23066 cals. per volt. 
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Fig. 3. Ellingham diagram for chlorides based on the 
decomposition potential measurements of Camsr and 
Devoro. 


Fig. 2. Ellingham diagram for chlorides based on current 
potential measurements of NeuMANN and RicurTer. 


Table 3. Calculated approximate free energy of formation of various chlorides 





‘ I 
AS AG@™ 
AH™* ane , . _ 
o- 


Substance (< Q)) ] y ; . : S¥c1 Ss Mu l 
» FP rr é n A . 298 T 208 
(keala) . v.U. : (AH TAS™) 
(keals) (keals) 


LiCl 97-65 7 26-66 0-46 91-9 
NaCl 98-33 17) 2. 26-66 21-86 918 
KC! 104-36 19-35") Bet 26-66 22. 97-7 
MgCl, 53-3 21-4) 76 53-32 39-68 70-7 
CaCl, 190-6 27-2) ot 53-32 36-02 x09 
SrCl, 197-87 27-5) 3° 53-32 30-12 93-1 
BaCl, 205-28 28-6 5 53-32 39-8: 96-7 
BeCl, 11246 We) 2-: 53-32 36-6 50-8 
CeCl, 250-8) 42°) 79-98 51-7! 81-4 
LaCl, 263-7 we 13-7 79-98 516s 82-8 
ThCl, 284-5) 54) 13-6 106-64 iti: 66-2 
AICI, ;,) 167-9) 22-4(°) 6-73 79-98 Me 49-6 
ZC, 268-9 a5 9-3 106-64 61-9 





* Unless otherwise stated the figures have been taken from Brcnowsk1 and Rossin1 (op. cit.). 

t Unless otherwise stated the figures have been taken from WENNER (op. cit.). 

‘ Bower and Houmann [69]. 

(») B.H. Zowm and J. E. Maver [70]. 

‘ Kewiey [71]. 

‘) K.K. Kerrey and G. E. Moore (72). 

‘©’ Estimated from a plot of entropy against log mol. wt. for equi-valent halides based on data given by Wenner, op. cit. 
and subsequent literature. 

() Rors and Biicuner [73]. 

() Metals Reference Book [64]. 

G) Estimated from Larmers formula [63]. 

(*) K. K. Kecury (74). 

() 8.8. Topp [62]. 

‘ L. Eyvrine and E. F. Westrum [76). 
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Table 4. Calculated approximate temperature variation of free energy of formation of various chlorides 





| Phase transformations AG|T equation: AG = A+ B(T—C 

AH™ AS™ , ! ne 

Substance Table 3 Table 3 
(keals) (B.U.) 


AStrans 
(=L/T) 
(BE.U.) 


Temp. Nature of 


Range of T. . l 
((K) transform. A 1000 B C 


(-K) n 


97-65 298- 459 97-65 19-46 
Li melts 1-66 459— 887 88-72 21-12 
887 LiCl melts 5-63 887-1300 79-68 15-49 


98-33 298— 371 98-33 21-86 
371 Na melts 1-70 371-1074 9-22 23-56 

1074 NaCl melts 6-73 1074-1151 73-66 16-83 

1151 Na boils 20-2°¢) 1151-1300 72-36 37-03 

298— 337 104-36 22-51 

337 K melts 1-92 337-1031 96-78 24-43 

1031 K._ boils 204) 1031-1063 79-82 44-43 

1063 KC] melts 6-02 1063-1300 78-4 38-41 

39-68 208— 924 153-3 39-68 

924 Mg melts 1-95 924— 986 116-64 41-63 

986 MgCl, melts 10-44" 986-1300 114-06 31-19 

190-6 298-1047 190-6 36-02 
1047 CaCl, melts 642 1047-1124 15236 29-60 

1124 Ca melts 184° 1124-1300 150-61 31-44 


197-87 298-1044 197-87 39-12 
1044 Sr melts 1-37 1044-1145 157-03 40-49 

1145 SrCl, melts 6-5 1145-1300 152-04 33-99 

205-28 208— 977 205-28 39-82 
977 Ba melts 1-51 977-1233 166-38 41-33 

1233 BaCl, melts 1233-1300 155-8 J4-83 

298-1088 259-8 51-78 

1088 Ce melts 1088-1121 203-46 54-28 

1121 CeCl, melts ef 1121-1300 201-67 47-78 





* See note * to Table 2. 

(») G. E. Moore [75). 

(9 Metals Reference Book [64]. 

(4) Estimated from known values for similar elements and compounds. 
(©) S. GuassTone, op. cit., p. 449. 


and 4, the calculations being carried out in the approx- In comparing the figures of NEUMANN and Ricu- 
imate manner already used for fluorides. TER, CamBt and Devoro, and the present rough cal- 


Probable electropotential order for the chlorides culations, several features are seen to be common to 


(a) Alkali and alkaline earth metals at least two of the three sets of results. 


In view of the various phase changes of the salts and 1. The most electropositive metal is barium. 
their constituent elements which occur within the 


goat ; 2. Potassium follows next at the lower tem- 
temperature range of practical interest and the in- ; me 
: - ‘ peratures but becomes relatively less electropositive 
stances of crossing of the decomposition potential ; : 
; Mie ; as the temperature increases. 
temperature curves found by various workers, the 


chloride electropotential order may be expected to 3. The KCl and LiCl curves cross between 1100 
show considerable changes with temperature. and 1250° K. 
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Table 5. Calculated approximate temperature variation of free energy of formation of some oxides 





AsS™ 
(=Syo ! 


Sy—So) ,, 


oo +?oe 
Ss MO s M 


728 


AH™ 
(keals) 


For No. of gm. atoms in formula 


weight (E.U.) (°K) 


Li,O 1423 906t 134 669 459 


K,O 337 


86-2 20-6* 304 38-0 1031 


BeO 135-0 3-37 64-0 


BaQO 133-0 16-8 63-1 977 


La,W, 457-0 30* TO-95 724 1158 


CeO, 233-4 17-7 49-03 45-13 55-0 1088 


Tho, 293-0 19-6 = 13-6TF 49-08 43-03 TOO 


SrO 140-8 13-0 13-3 24-52 24-82 66-7 


Phase transformations 


(keals) Temp. Nature of 


AG/T equation: AG = A + B(T — C) 


1¢ (1300 


Range 
of T. 


! St rans 


transform. (E.U.) L000 B C on 


(keals) 


142°3 
120-1 


28-86 0 
30-52 459 


208 
459 


459 


Li melts 1300 


86-2 
74-64 
49-49 


34-32 0 
36-24 337 
56-24 1031 


208 
337 
1031 


337 
1031 
1300 


K melts 
K. boils 
135-0 28:43 


208-1300 


133-0 
110-705 


22-82 


24-33 


208 


977 


977 


1300 


Ba melts 


FO-O5 0 
73-45 1158 


457-0) 
374-84 


208-1158 


1300 


La melts 

1158 
45-13 0 
47-63 1088 


233-4 
184-3 


208 LOSS 
1088-1300 


Ce melts 


208-1300 293-0 43-03 0 


1300 140-8 24-82 0 





Notes. Except where otherwise stated values of 1H have been taken from Brcnowsk1 and Rossini (op. cit.), 


WENNER (op. cit.) and other data from the previous tables. 


entropies from 


* Estimated from plot of entropy against log mol. wt. for oxides of metals of similar valency. 


** Estimated from known values for similar elements. 
+ Jounstone and Baver [77]. 
tt Metals Reference Book [64]. 


4. The NaCl and CaCl, curves lie close together 
and possibly intersect*. 

5. Sodium and calcium are the least electropositive 
of the metals considered. 

If the slopes of the free-energy temperature curves 
are considered, it will be noticed that Neumann and 
RicuTer found the chlorides to fall into well defined 
groups, the members of which had sensibly the same 
slope: the alkaline earth metal chloride curves had 
a comparatively gentle slope, LiCl was steeper, and 
the alkali metal chlorides had a uniformly steep slope**. 
Such a grouping of the slopes was not however found 
by Campr and Devoro, nor does it appear probable 
from the present rough calculations. 

None of the experimental determinations of dE/dT 
indicate the changes of slope which would be expected 
to occur at the temperatures corresponding to phase 
changes of the chlorides or their constitutent elements 

* An intersection of these curves between 600 and 800°C 
would accord with the observations of DanneeL, STockEM 
and v. KUEGELGEN [28] in electrolysing NaCl-CaCl, mixtures. 


** NeumMaNN and Ricurer include figures for RbCl and 
CsCl. 


and which in the case of KCl appear to be of appreci- 
able magnitude. 

With regard to the actual values of the slopes of the 
free-energy temperature curves, the calculated slopes, 
which should presumably be at least of the right order, 
agree quite closely with the lowest experimental 


slopes (Cambri and Devoro) in both the cases which 


have been fairly thoroughly investigated, namely, NaCl 
and KCl, and where the data available for the rough 
calculations are probably the most reliablet. In addi- 
tion the calculated free energies of formation agree 
quite closely with the measurements of Cams and 
Devoro (Table 6). 
AGT curves for NaCl and KCI can probably therefore 


be obtained by plotting the mean values of AG and 

(AG) 
eT 

The general form of the curves for the other chlor- 


ides considered can then be roughly indicated by 


A close approximation to the 


given in the table (Fig. 4). 


t In these two cases, moreover, complications due to 
incomplete dehydration and sub-chloride formation do not 
arise. (Krowt [29] states that in the electrolysis of BaCl, 
sub-chloride only is produced.) 
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0) for NaCl and KCl (keal) 





Substance 
Temp. 


(°K) 


1100 
1200 


1000 


(AG) 
c a 


c 


Temp. 


( K) Mean 


Mean 


1150 


> 1050 





Campi and Devoro (22); EF 


relating them with those of NaCl and KC! in accor- 
dance with the broad conclusions already reached and 
giving them as far as possible the calculated gradients, 
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Fig. 4. Ellingham diagram for chlorides based on all data 
discussed. 


consistent with their passing through the calculated 
values at 208° K [30). 


(b) Magnesium and the rare earth metals 

Various determinations have been made of the de- 
composition potential of magnesium chloride [4], {10}, 
[31 }-[34]. Without attempting a critical discussion 


of these figures, which range from 2-50 to 2-59 vy at 
the melting point with temperature coefficients vary- 


Calculated figure from Table 4. 


ing from 0-00057 to 0-00000 vy deg.. the most satis- 
factory measurement appears to be that of Conen, 
not only because of the care taken in the deter- 
mination but also because his value of 13 cals deg. 
C for dE dT, which is the lowest, agrees fairly well 
with caleulation (15-5 cals deg.). Lorenz and VELDEs 
value of dE dT is almost identical with the calculated 
figure, but their value of EF is much the lowest and 
likely to be low from consideration of the experi- 
mental procedure*. Conens values for both EF and 
dE dT are therefore used for plotting on Fig. 4, the 
value for 1G°*** being taken from Table 3. 

No reliable determinations of the decomposition 
potential of the rare earth metal chlorides appear to be 
available. NeuMANN and Ricurers figures are cert- 
ainly low in absolute value, ranging at 800°C from 
1-3 v for PrCl, to 1-9 v for CeCl,. They are also very 
probably low relatively to MgCl, which is given as 
2-2 v instead of 2-5 v; for it is well known in the 
industry that cerium will reduce MgCl, to a consider- 
able extent, MgCl, being in fact usually excluded from 
fluxes for alloys containing rare earth metals [35 }- [37 |. 
The value of dE dT found by Neumann and RicuTer 
(0-002 vy °) moreover, corresponds to an entropy of 
formation for CeCl, of 138 £.U'. compared with the 
calculated figure of 52. Drosssacn has studied the 
electrolysis of the double chlorides of cerium and 
neodymium with potassium, and the deposition 
potentials observed (Table 7), although probably not 
simply related to the reversible potentials, appear to 
indicate (a) that the reversible potentials are similar 
to but higher than that of MgCl, thus according quali- 
tatively with calculation**, and (b) that (in the chloride 

* The chlorine stream was interrupted during the measure- 
ment and the voltage observed remained constant for only 
a very short interval. 

** Tho caleulated free energy difference corresponds to 
nearly half a volt. 
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Table 7. Deposition potentials observed by Drosspacn 





Reversible 


decomposition 


Temp. Deposition 
Sale (C) votential (v) potential (v) 
Mer 
’ of single chloride 
where known 
2KCIl MgCl, 740) 2-71-2-85 2-54 (740° ©) 
247 (860°C) 
2KCl CeCl, S60) 2-96-3-10 
2KCl NdCl, 860 2-75-3-00 





series) cerium is a little more electropositive than neo- 
dymium, as was found by Neumann and RicurTer. 
Experiments on the reduction of NdCl, and GdCl, by 
magnesium have been made by TromBe and Maun [38] 
who reproduce NeuMANN and RIcHTERs measure- 
ments without comment as to the probable accuracy 
of the decomposition potentials and their temperature 
coefficients. Trompe and Maun obtained percentage 
reductions of the order of 50%, indicating either that 
magnesium and the rare earth metals do in fact differ 
little in electropotential in their chlorides, or, more 
probably, that the activity of cerium in magnesium 
is low. This would be expected in view of compound 
formation. The small difference in maximum per- 
centage reduction obtained with NdCl, and GdCl, 
(42 and 52% respectively) is such as to cast doubt on 
the large difference in decomposition potential of 
0-65 v between CeCl, and NdCl, found by Neumann 
and Ricurer. A difference of this magnitude should 
render fractionation of the rare earth metals by 
electrolysis or chemical exchange a simple process. But 
in fact no appreciable enrichment in any constituent 
can be obtained by (a) reducing mischmetal* chloride 
or fluoride with magnesium (b) ‘‘fluorinating’’ misch- 
metal chloride with KF or (c) fluxing magnesium- 
mischmetal alloys with fluxes containing MgC\,. 
The 


chlorides and for MgCl, should lie very close together 


1G T curves for the various rare earth metal 


and slightly below** (i.¢. representing higher numerical 
free energy) that of MgCl, (Fig. 4). 
NEUMANN 


In accordance 


with and RicuTers measurements, and 


Trompe and Manns observations that the percentage 


* “Mischmetal”’ refers to a mixture of the cerium group 


of rare earth metals roughly in the proportions Ce 50%, 
La 25%, Nd 18%, Pr 7%. 

** Even if Raou ts law is assumed, Trompe and MAHNs 
percentage reduction figures of about 50% with equal weights 
of magnesium and “REC,” (i.e. a 10:1 molecular ratio) 
would correspond to values of K (the equilibrium constant) 
less than unity for the reaction 3Mg + 2RECI, — 3MgCl, 4 
2RE, 
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reduction increases with temperature, the slopes of 
the AG T curves for the RE metals are probably 
somewhat steeper than that of the MgCl, curve, the 
contrary indication of the calculation being here of 
little value in the absence of measured values for the 
entropy of CeCl, and the latent heats of fusion of 
cerium and CeCl,. 

have been 


Magnesium reduction 


carried out on ThCl, by Fox [39], but the low alloying 


experiments 


efficiencies obtained do not necessarily indicate that 
the salt is only slightly reducible since the sample used 
may not have been completely anhydrous. The rough 
calculations indicate that reduction of ThCl, by mag- 
nesium should be substantially complete and this has 
been confirmed experimentally. It is therefore evident 
that in the chloride electropotential series thorium lies 
somewhat below magnesiumt. Zirconium must lie 
well below magnesium since reduction of ZrCl, by 
magnesium is complete and the decomposition poten- 
tial of the salt appears to be very low [40]. 


Il. THe Oxipe ELecrroporentiaAL SERIES 
A recent critical compilation of free energy data of 
oxides has been made by RicHarpson and JEFFES [42]. 






























The results of rough thermodynamical calculations are 
given in Table 5 for such of the oxides not considered 
by RicHarpDson and Jerres as are of present interest. 
These data are combined in Fig. 5. Since preparing 
this paper, calculated data have been published for 
BaO and ZrO, by Ossporn [43]. His results for BaO 
agree very closely with Table 7. Those for ZrO, have 
been added to Fig. 5, but it must be remembered that 
they are not of much significance where reduction is 
concerned in view of the readiness with which oxygen 
dissolves in zirconium. 


IV. Nore oN APPLICATION OF THE FREE ENERGY 
Data TO MaGNesium TECHNOLOGY 
The free energy data collected in Figs. 1, 4 and 5 should 
enable very useful indications to be obtained immedi- 
ately of the practicability of proposed reactions tf. 
Strict application of the data to a particular case, 
would however require additional knowledge, such 
as the free energies of formation of compounds 
+ This accords with Haissinskys claim [41] that thorium 
is a true analogue of zirconium and resembles it in electro- 
chemical! properties ; but in the fluoride series thorium certainly 
resembles the rare earth metals, ThF, being very incompletely 
reduced by magnesium. 


++ For examples of the application of such data see 
references [42], [44], [45]. 
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formed in the metal or flux phases and the relevant 
activity coefficients*. While the experimental deter- 
mination of such data might yield results of great 
interest from a fundamental point of view, the practical 
value of the work would be small in proportion to the 


Je mperature ——= 


700 _ _ 900 





= 
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Temperoture ———e C 


5. Ellingham diagram for oxides. (The curves for 
Al, Mg and Ca are due to RicHarpson and Jerres 
that for Zr to OsBorn. 


thermochemical calculations.) 


Others are based on rough 


labour involved. Of greater practical value, however, 
would be a knowledge of the concentrations of sodium, 
potassium, etc. in magnesium which are in equilibrium 
with constant molecular proportions of their chlorides. 
Such information would enable one to draw up a 


“Chloride Chemical Displacement Series” with re- 
* It is often possible to infer qualitatively the effect of the 
deviations from RaouLtT: law to be expected in particular cases. 
For example sodium shows limited miscibility with magnesium 
and MgCl, forms a compounds with NaCl. These deviations 
are likely to offset each other, so that the sodium content of 
magnesium made electrolytically from a bath containing 
NaC! is likely to be of the order indicated by Raovu ts law. 
On the other hand calcium forms a compound with magnesium 
and the two chlorides from a eutectic system, so that the 
calcium content of electrolytic magnesium is likely to be 
onsiderably greater than Raovutt: law indicates. These 
onclusions are borne out by published data on the sodium 
d calcium contents of electrolytic magnesium [3] and the 
omposition of cell electrolyte [50]. 
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spect to magnesium which would be of immediate 
practical application. The order of elements in this 
chemical displacement series might differ considerably 
from that in the electropotential series in the case of 
elements with chlorides of similar stability. Thus the 
free energy curves for NaCl and LiCl are seen from 
Fig. 4 to lie close together and, in the absence of 
activity data, it would be thought that NaCl would 
tend to flux lithium from magnesium**. Actually the 
reverse tendency is found to occur, LiCl being quite 
effective in removing sodium [46]. This could be fore- 
seen at once from a chemical displacement series. 
The exchange experiments needed to establish the 
latter series would require precautions to ensure the 
achievement of equilibrium; and in determining a si- 
milar series for fluorides it might be difficult to avoid 
complicating the issue by using a suitable solvent for 
the fluoride f. 
be used in industry, so that the thermal effects of 


Such solvents would however have to 


possible double halide formation would not invalidate 
the results from a practical point of view tt. From the 
experimental results obtained, rough estimates of the 
amounts of impurities in metal treated with any 
particular flux composition could be made. 


DousLeE Decomposirion Reactions IN MOLTEN 
Sat MIXTURES 

Since in fused salt mixtures the constituents appear 
to be completely ionised [47], the contrary evidence 
of Sackur [48] from freezing point determinations 
resting on an erroneous basis [49], it may be assumed 
that the compounds present at equilibrium will depend 
on the activity solubility products of the various 
The insolubility of MgF, in 
chloride mixtures other than those very rich in MgCl, 


possible compounds. 


has been mentioned elsewhere [50], and in many ex- 
periments with fluxes and alloying salts it has been 
observed that, wherever MgF, could form from the 


** Recent determinations of the decomposition potential 
of LiCl by Grorue and SAVELSBERG [34] actually give lower 
values than that here accepted for NaCl. The temperature 
coefficient found is however twice the probable value, and, 
since the decomposition potential found for KCl is also con- 
siderable lower and its temperature coefficient higher than the 
probable values, these measurements evidently do not rep- 
resent true reversible values. 

+ The alkaline earth metal fluorides are solid at normal 
alloying temperatures (700-850° C) and would not form with 
the MgF, produced on reaction mixtures melting within this 
temperature range. 

++ Some work has in fact been carried out by the author 
along these lines, but the results are not yet complete. 
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constituents present, it did so and was precipitated out, 
Rare earth metal fluorides are also insoluble in most 
chloride mixtures, and, in the absence of MgCl,, any 
soluble fluoride will generally precipitate as “REF,” if 
rare earth metal chlorides are present. The solubility 
of other fluorides in chloride mixtures depends very 
much on the composition of the latter. In CaCl, 

NaCl— KCl mixtures, for example, BaF, and KF are 
relatively insoluble and NaF and CaF, relatively 
soluble. The stability at lower temperatures of KF 
in certain salt mixtures containing CaCl, and NaCl 
51) is an example of the effect of solubility governing 
the constituents present in molten salts, since on ther- 
mochemical grounds CaF, would be expected, and even 
NaF should be formed in preference to KF. The pre- 
ferential precipitation of fluorine as MgF, and “REF,” 
is however to be expected on thermochemical as well 
as solubility considerations, and it is therefore of inter- 
est to consider what metathetic reactions would be 
expected to oceur from thermochemical calculations 
alone, assuming complete solution of all constituents, 
and then to see how these results compare with 


observation in a few test cases. 


CONVERSION OF FLUORIDE TO CHLORIDE 


Consider the conversion of fluoride to chloride re- 
ZCl MF. 


The change in free energy per mol (4@) may be 


presented by the equation: MCI — ZF 


regarded as being composed of two terms 


( Gy) Miyy) and ( Ivy yy, ) 


Since for any given fluoride, ZF, the first term is 
constant, the ability of chlorides to convert fluorides to 
chlorides will depend on the value of the second term. 
The second term will however always be of the same 
sign whether the reaction tends to go from left toright 
or vice versa. To obtain a rough measure of the relative 
“chlorinating powers” of the various chlorides, a stand- 
ard fluoride could be chosen, preferably one whose 
chloride did not possess a high heat of solution so that 
the relative positions of the metal in the chloride 
and fluoride electropotential series would be similar. 
But it is probably more satisfactory to express the 
chlorinating powers relative to an “‘ideal’’* fluoride by 
evaluating the term (A4G,y—AGy) and adding 
this to the second term (AG yy, 
the chlorinating powers of various chlorides calculated 


IG, \)- Table Ss gives 


* This would possess zero heat of solution as would the 
corresponding chloride. 
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in this way from the 14 values at room temperature ** 
given in Tables 1-4. 
Table 8. chloride 


Conversion of fluoride to approximate 


chlorinating powers of various chlorides 





Chlorinating 
Chloride power + 


Col. 1 + Col. 2 


l 
1G", —AG™..)\ AGry 1p: 
_ AO ur 4c) 4Scr,, 49 Pr, 


BeCl, 20-7* 
MgCl, 
AICI, 19-6 
CeCl, 15-3* 
CaCl, t 15-1 
LiCl 13-2 
BaCl, > + 33- 7-4 
NaCl + 33+ 3-2 
KCl 20) 33°! » 40 


21-0 





+ i.e. Free energy increase (keal) per gram equivalent 


accompanying conversion of an ‘ideal’ fluoride to an “ideal ’’- 

chloride. (n = valency of cation.) 

* The correct figure will be numerically less, the value of 

+ 2 HP in the 
9 


absence of information on the value of the heat of solution. 


1Hir, having been here taken as THRs 


The same applies to JH@%p.. 


The arrangement of the common Group IT chlorides 
in the table, namely, Mg—Ca—Ba in order of de- 
creasing chlorinating power, agrees with that observed 
in the chlorination of ZrF, [15]. The table also shows 
clearly why KCl will not chlorinate ZrF,; and why in 
and K°Ca’ 


the stable fluoride is in each case CaF, [52}. 


the reciprocal systems Na‘’Ca’* Cl F 
Ci-F-, 
The correctness of the relative positions of BaCl, and 
LiCl indicated in the table has been checked experi- 
mentally, using a previously described method [50] 
to avoid possible interchange of radicals during analysis 
of the fused salt mixture obtained. 


Conversion of oxide to chloride 


Approximate date for conversion of oxides to chlorides 
are given in Table 9, the chlorinating power being here 
expressed as the free energy increase per gram equi- 
valent accompanying conversion of BaO to BaCl, 
by the chloride in question. The barium compounds 
are taken as reference standard, since BaCl, has the 
smallest heat of solution of any of the chlorides in the 
table, with the exception of NaCl, in which case the 
chloride is obviously relatively much more stable than 

** The calculated order of chlorinating power at 1000°K 


is unaltered so far as concerns the chlorides for which high 
temperature data have been estimated. 
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Table 9. Conversion of oxide to chloride—approximate relative 
chlorinating powers of various chlorides 





Approximate 
relative 
(46% AGYa) : ( IGRaci.—4Ghao) chlorinating 
- power * 
Col. 2 +- Col. 3 


Chloride 


46-8 
44-9 
32-1 


33-6 
33-6 
33-6 
33-6 27-0 
33-6 15-8 
33-6 9-6 
33-6 8-9 
33-6 0 

33-6 13-7 


KCl 33-6 ~ 26-1 





* i.e. Free energy increase (keal) per gram equivalent 
accompanying conversion of BaO to BaCl,. (nm = valency of 


cation.) 


the oxide. The order of ability to chlorinate oxides found 
by calculation for the alkali and alkaline earth metal 
chlorides is probably correct, since experiments on the 
addition of BaO to melts of MgCl,—CaCl,—NaCl— KC] 


Table 10. Comparative values of free energy of formation per gram 

: { (71000 

equivalent at 10007 K jestimated by various workers 
n 





Fluoride Chloride 


Di ife rences Dit fe rencea 


B-—E K—E* B—E K—E* V—E 


122-2 +2: 0 7+ +21 
oo +18 LOS 5-6 -O7 
109-4 5 1 OS ‘ 5 ‘ 1-1 
110-3) «4+-2- LOS 1-1 
1240 +04 L110 75-6 2-6 
124-0 N 79-4 OS 
123-6 +0. -O5 S28 +01 





* Differences are here given only to the nearest half kilo- 
‘al. in view of the difficulty of reading the diagram given by 
KELLOGG to a greater accuracy. (2 Estimate from Table 2 
w Fig.4; B — Estimate of Brewer, BromMLey, GriLies and 
Lorcren; A — KeLioaes estimate; V =VILLAs estimate; O 
OSBORNS estimate; N = No estimate made.) 


‘lectrolyte showed that magnesium was first preci- 
pitated and only when the MgCl, had been reduced 
o below 1% was calcium precipitated. No Na,O or 


r 


20 were precipitated. These results would also be 
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expected to occur on grounds of solubility since MgO 
is apparently insoluble in all mixtures of alkali and 
alkaline earth metal chlorides [53], CaO is soluble 
in CaCl,, while SrO and BaO appear to show wider 
solubility [48]. 


The author is indebted to Dr. O. 
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